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METHOD AND SYSTEM FOR EVALUATING 

CARDIAC ISCHEMIA 



10 



15 



20 



Joseph M. Starobin and Yuri B. Chemyak 



• Related Applications 

This application is a continuation-in-part of copending application serial 

number 09/603,286, filed June 26, 2000, the disclosure of which is incorporated by 

> ■ 

reference herein in its entirety. 



Field of the Invention 

The present invention relates to non-invasive high-resolution diagnostics of 

■ 

cardiac ischemia based on processing of body-surface electrocardiogram (ECG) data. 
The invention's quantitative method of assessment of cardiac ischemia may 
simultaneously indicate both cardiac health itself and cardiovascular system health in 
general. 



Background of the Invention 

Heart attacks and other ischemic events of the heart are among the leading 

r 

25 causes of death and disability in the United States. In general, the susceptibility of a 

4 

particular patient to heart attack or the like can be assessed by examining- the heart for 

i 

evidence of ischemia (insufficient blood flow to the heart tissue itself resulting in an 
insufficient oxygen supply) during periods of elevated heart activity. Of course, it is 
highly desirable that the measuring technique be sufficiently benign to be carried out 
30 without undue stress to the heart (the condition of which might not yet be known) and 

« ■ 

without undue discomfort to 1 the patient. \ 

Hie cardiovascular system responds to changes in physiological stress by 
adjusting the heart rate, which, adjustments can be evaluated by measuring the surface 
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■ 

ECG R-R intervals. The time intervals between consecutive R waves indicate the 
intervals between the consecutive heartbeats (RR intervals). This adjustment normally 
occurs along with corresponding changes in the duration of the ECG QT intervals, 
which characterize the duration of electrical excitation of cardiac muscle and 
5 represent the action potential duration averaged over a certain volume of cardiac 
muscle (Figure 1). Generally speaking, an average action potential duration measured 
as the QT interval at each ECG lead may be considered as an indicator of cardiac 

4 

systolic activity varying in time. 

Recent advances in computer technology have led to improvements in 
10 automatic analyzing of heart rate and QT interval variability. It is well known now 
that the QT interval's variability (dispersion) observations performed separately or in 
combination with heart rate (or RR-interval) variability analysis provides an effective 

- 

■ » * 

tool • for the assessment of individual susceptibility to cardiac arrhythmias 

4 

(RSurawicz, J. Cardiovasc. Electrophysiol, 1996, 7, 777-784). Applications of 
15 different types of QT and some other interval variability, to susceptibility to cardiac 
arrhythmias are described in U.S. Patents by Chamoun No.5,020,540, 1991; Wang 

* * 

No. 4,870,974, 1989; KroU et al. No.5,1 17,834, 1992; Henkin et al. No. 5,323,783, 
1994; Xue et al. No.5,792,065, 1998; Lander No.5,827,195, 1998; Lander et al. 
No.5,891, 047, 1999; Hojumetal. No.5,95 1,484, 1999). 

# 

20 It was recently found that cardiac electrical instability can be also predicted by 

linking the QT - dispersion observations with the ECG T-wave alternation analysis 
(Vender et al., U.S. Patents No.5,560,370; 5,842,997; 5,921,940). This approach is 
somewhat useful in identifying and managing individuals at risk for sudden cardiac 
death. The authors report that QT interval dispersion is linked with risk for 

25 arrhythmias in patients with long QT syndrome. However, QT interval dispersion 
alone, without simultaneous measurement of T - wave alternation, is said to be a less 
accurate predictor of cardiac electrical instability (U.S. Pat. 5,560,370 at column 6, 

■ 

lines 4-15). 

» » 

Another application of the QT interval dispersion analysis for prediction of 
30 sudden cardiac death is developed by J. Sarma (U.S. Patent No. 5,419,338). He 
describes a method of an autonomic nervous system testing that is designed to 
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* 

evaluate the imbalances between both parasympathetic and sympathetic controls on 
the heart and, thus, to indicate a predisposition for sudden cardiac death. 

The same author suggested that an' autonomic nervous system testing 

» 

procedure might be designed on the basis of the QT hysteresis (J.Sarma et al., PACE 
5 10, 485-491 (1988)). Hysteresis between exercise and. recovery was observed, and 
was attributed to sympatho-adrenal activity in the early post-exercise period. Such an 

* ■ 

activity was revealed in the course of QT interval adaptation to changes in the RR 
interval during exercise with rapid variation of the load. 

The influence of sympatho-adrenal activity and the sharp dependence of this 
10 hysteresis on the time course of abrupt QT interval adaptation to rapid changes in the 

* * * • 

' RR interval dynamics radically overshadows the method's susceptibility to the real 
ischemic-Iike changes of cardiac muscle electrical parameters and cardiac electrical 

» 

conduction. Therefore, this type of hysteresis phenomenon would not be useful in 

* 

assessing the health, of the cardiac muscle itself, or in assessing cardiac ischemia. 
15 A similar . sympatho-adrenal imbalance type hysteresis phenomenon was 

+ 

observed by A. Krahn et al. {Circulation 96, 1551-1556 (1997)(see Figure 2 therein)). 
The authors state that this type of QT interval hysteresis may be a marker for long-QT 
syndrome. However, long-QT syndrome hysteresis is a. reflection of a genetic defect 
of intracardiac ion channels associated with exercise or stress-induced syncope or 
20 sudden death. Therefore, similar to the example described above, although due to two 

• * 

different reasons, it also does not involve a measure of cardiac ischemia or cardiac 
muscle ischemic health. 

A conventional non-invasive method of assessing coronary artery diseases 

* • 

associated with cardiac ischemia is based on the observation of morphological 

« 

25 changes in a surface electrocardiogram during physiological exercise (stress test). A 
change of the ECG morphology, such as an inversion of the T-wave, is known to be a 
qualitative indication of ischemia. The dynamics of the ECG ST- segments are 
continuously monitored while the shape and slope, as well as ST-segment elevation or 
depression, measured relative to an average base line, are altering in response . to . 

30 exercise load. A comparison of any of these changes with average values of . 
monitored ST segment data provides an indication of insufficient coronary blood 
circulation and developing ischemia. Despite a broad clinical acceptance and the 
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availability of computerized Holter monitor-like devices for automatic ST segment 
data processing, the diagnostic value of this method is limited due to its low 
sensitivity and low resolution. Since the approach is specifically reliable primarily for 

» 

ischemic events associated with relatively high coronary artery occlusion, its 
5 widespread use often results in false positives, which in turn may lead to unnecessary 

and more expensive, invasive cardiac catheterization. 

. ■ 

Relatively low sensitivity and low resolution, which are fundamental 
■ disadvantages of the conventional ST-segment depression method, are inherent in 
such method's being based on measuring an amplitude of a body surface ECG signal, 

10 which signal by itself does not accurately reflect changes in an individual cardiac 
cell's electrical parameters normally changing during an ischemic . cardiac event. A 
body surface ECG signal is a composite determined by action potentials aroused from 
discharge of hundred of thousands of individual excitable cardiac cells. When 
electrical activity of excitable cells slightly and locally alters during the development 

15 of exercise-induced local ischemia, its electrical image in the ECG signal on the body 
surface is significantly overshadowed by the aggregate signal from the rest of the 
heart. Therefore, regardless of physiological conditions, such as stress or exercise, 
conventional body surface ECG data processing is characterized by a relatively high 
threshold (lower sensitivity) of detectable ischemic morphological changes in the 

9 

20 ECG signal. An accurate and faultless discrimination of such changes is still a 
challenging signal processing problem. 

Accordingly, ah object of the present invention is to provide a non-invasive 
technique for detecting and measuring cardiac ischemia in a patient. 

• ■ 

Another object of the invention is to provide a technique for detecting and 
25 measuring cardiac ischemia, which technique is not unduly uncomfortable or stressful 
for the patient. 

Another object of the invention is to provide a technique for detecting and 
measuring cardiac ischemia, which technique may be implemented with relatively 
simple equipment. 

■ * * 

30 Still another object of the invention is to provide a technique for detecting and 

» 

measuring cardiac ischemia, which technique is sensitive to low levels of such 

■ 

ischemia. 
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Summary of the Invention 

The present invention overcomes the deficiencies in the conventional ST- 

9 * 

segment analysis. Although still based on the processing of a body surface ECG 
5 signal, it nevertheless provides a highly sensitive and high resolution method for 
distinguishing changes in cardiac electrical conduction associated with developing 
cardiac ischemia. In addition to the significant cardiac ischemic changes detectable by 

* ■ ■ 

the conventional method, the present invention allows one to. determine much smaller 
ischemia-induced conditions and alterations in cardiac electrical conduction. Thus, 

* > 

10 unlike a conventional ST-segment depression ischemic analysis, the method of the 
present invention opens up opportunities to detect low-level cardiac ischemia 
(undetectable via the regular ST-segment method) and also to resolve and monitor 
small variations of cardiac ischemia. In particular, individuals who would be 
considered of the same level of cardiac and cardiovascular health according to a 

9 

• m 

15 conventional ECG evaluation (an ST-depression method), will have different 
measurements if compared according to the method of the present invention, and the 
cardiac and cardiovascular health of an individual can be quantitatively evaluated, 

■ * * 

. compared and monitored by repeated applications of the method of the present 
invention. • 

20 • The present invention is based in part on the discovery that, under certain 

physiological conditions, QT- and RR- interval data sets may be interpreted as 
representing composite dispersion-restitution curves, which characterize the basic- 
dynamic properties of the medium (in this case, cardiac muscle). Indeed, if rapid QT 
interval adaptation facilitated by sympatho-adrenal activity occurs much faster than 

■ ■ ♦ 

25 gradual heart rate changes following slow alteration of external physiological 
conditions, then the QT interval may be considered primarily as a function of a heart 
rate and/or a preceding cardiac cycle length and does not substantially depend on 
time-dependent sympatho-adrenal transients. In such a case a: particular QT-.and KR- 
interval data set determines a time-independent, dispersion-like, quasi-stationary 

30 curve which does not substantially depend on rapid adaptational transients and 
depends primarily on medium electrical parameters. 



WO 02/00113 PCT/US01/20306 

-6- 

» 

Based on this discovery, the present invention provides a highly sensitive and 
high resolution method of assessing cardiac ischemia. This method allows one to 
detect comparatively small alterations of cardiac muscle electrical excitation 
properties that develop during even a moderate ischemic condition. For example, 

* m 

5 consider a gradual heart rate adjustment in a particular human subject in response to 
slow (quasi-stationary), there-and-back changes of external physiological conditions. 
Ideally, when a cardiac muscle is supplied by a sufficient amount of oxygen during 
both gradually increasing and gradually decreasing heart rate stages, the 
corresponding, there-and-back, quasi-stationary QT versus RR interval curves which 
10 result should be virtually identical. However, if ischemia exists, even if only to a very 
. . minor extent, there will be alterations of cardiac muscle repolarization and excitation 
properties for the human subject with the result that one observes a specific quasi^ 

* * 

stationary QT-RR hysteresis loop. Unlike non-stationary loops (J. Saima et ah, supra 
(1987); A. Krahn et al., supra (1997)), the quasi-stationary hystereses of the present 

15 invention do not vary substantially versus the course of sympathoadrenal QT- 
k interval adjustment. The domains and shapes of these loops are not significantly 
affected by time-dependent transients rapidly decaying during a transition from one 
. , particular heart rate to another; instead, they depend primarily on ischemia-induced 
changes of medium parameters. The domain encompassed by such a quasi-stationary 

20 hysteresis loop and its shape represent a new quantitative characteristics that indicate 
cardiac muscle health itself and the health of the cardiovascular system in general. 
Moreover, any measure of the shape and/or domain enclosed in the hysteresis loop (a 
measure of a set as defined in the integral theory) possesses the property that any 

■ 

expansion of the domain results in an increase of the measure. Any such mathematical 
25 measure can be taken as the new characteristics of cardiac health mentioned above. 

4 * 

An arbitrary monotonic function of such a measure would' still represent the same 
measure in another, transformed scale. 

■ 

A first aspect of the present invention is a method of assessing cardiac 
ischemia in a subject to provide a measure of cardiovascular health in that subject. 
30 The method comprises the steps of: 

(a) collecting a first QT- and RR- interval data set from the subject 
* during a stage of gradually increasing heart rate; 
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(b) collecting a second QT- and RR- interval data set from the subject 
during a stage of gradually decreasing heart rate; ' 

(c) comparing said first QT- and RR- interval data set to the second 
QT- and RR- interval data set to determine the difference between the data 

5 sets; and 

(d) generating from the comparison of step (c) a measure of cardiac 
ischemia during exercise in said subject, wherein a greater difference between 
said first, and. second data sets indicates greater cardiac ischemia and lesser 

■ 

cardiovascular health in said subject. 
10 During the periods of gradually increasing and gradually decreasing heart rate the 
effect of the sympathetic, parasympathetic, and hormonal control on formation of the 
. hysteresis loop is sufficiently small, minimized or controlled so that the ischemic 
changes are readily detectable. This maintenance is achieved by effecting a gradual 
increase and gradual decrease in the heart rate, such as, for example, by controlling 
15 the heart rate through pharmacological intervention, by direct electrical stimulation of 
the heart, or by gradually increasing and gradually decreasing exercise loads. 

Accordingly, the foregoing method can be implemented in a variety of 
different ways. A particular embodiment comprises the steps of: 

(a) collecting a first QT- and RR- interval data set from said subject 

■ 

20 during a stage of gradually increasing exercise load and gradually increasing 

heart rate; 

♦ • 

(b) collecting a second QT-.and RR- interval data set from said subject 
during a stage of gradually decreasing exercise load and gradually decreasing 
heart rate; 

i ■ 

25 (c) comparing the first QT- and RR- interval data set to the second QT- 

and RR- interval data set to determine the difference between said data sets; 
and 

(d) generating from said comparison of step (c) a measure of cardiac 
ischemia during exercise in said subject, wherein a greater difference between 
30 said first and second data sets indicates greater cardiac ischemia and lesser 

cardiovascular health in said subject. 
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A second aspect of the present invention is a system for assessing cardiac 
ischemia in a subject to provide a measure of cardiovascular health in that subject. 
The system comprises: 

(a) an ECG recorder for collecting a first QT- and RR- interval data set 

• # 

5 from the subject during a stage of gradually increasing heart rate and 

r ■ 

collecting a second QT- and RR- interval data set from the subject during a 
stage of gradually decreasing heart rate; 

* 

(b) a computer program running in a computer or other suitable means 
for comparing said first QT- and RR- interval data set to the second QT- and 

* 

10 RR- interval data set to determine the difference between the data sets; and 

r « * 

t 

(c) a computer program running in a computer or other suitable means 
for generating from said determination of the difference between the data sets 
a measure of cardiac ischemia during exercise in said subject, wherein a 
greater difference between the first and second data sets indicates greater 

15 cardiac ischemia and lesser cardiovascular health in the subject. 

A further aspect of the present invention is a method of assessing cardiac 
ischemia in a subject to provide a measure of cardiac or cardiovascular health in that 
subject, the method comprising the steps, performed on a computer system, of: 

(a) providing a first QT- and RR- interval data set collected from the 
20 subject during a stage of gradually increasing heart rate; 

■ 

(b) providing a second QT- and RR- interval data set collected from 
the subj ect during a stage of gradually decreasing heart rate; 

(c) comparing the first QT- and RR- interval data set to the second QT- 
and RR- interval data set to determine, the difference between the data sets; 

25 and 

► 

(d) generating from the comparison of step (c) a measure of cardiac 
ischemia during stimulation in the subject, wherein a greater difference 
between the first and second data sets indicates greater cardiac ischemia and 

• * 

lesser cardiac or cardiovascular health in the subject. 
30 The first and second QT- and RR-interval data sets may be collected while 

minimizing the influence of rapid transients due to autonomic nervous system and 
hormonal control on the data sets. The first and second QT- and RR-interval data sets 
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are collected without an intervening rest stage. The generating step may be carried 
out by generating curves from each of the data sets, and/or the generating step may be 
carried out by comparing the shapes of the curves from data sets. In a particular 
embodiment, the generating step is carried out by determining a measure of the 
5 domain between the curves. In another particular embodiment, the generating step is 

■ 

carried out by both comparing the shapes of the curves from data sets and deterrmning 
a measure of the domain between the curves. The method may include the further 
step of displaying the curves. In one embodiment the comparing step may be carried 

• ■ 

out by: (i) filtering the first and second QT- and RRt interval data sets; (it) generating 
10 a smoothed hysteresis loop from the filtered first and second QT- and RR- interval 
data sets; and then (Hi) determining a measure of the domain inside the smoothed 

4 

hysteresis loop. In another embodiment, the comparing step may be carried out by: 
(i) filtering the first and second QT- and RR- interval data sets; (it) generating 
preliminary minimal values for the first and second QT- and RR- interval data sets; 
15 (Hi) correcting the preliminary minimal values; (iv) generating first and second 

■ 

preliminary smoothed curves from each of the filtered data sets; (v) correcting the 
preliminary smoothed curves; (vi) fitting the preliminary smoothed curves; (vii) 
generating a smoothed hysteresis loop from the first and second fitted smoothed 
curves; and then (viii) determining a measure of the domain inside the hysteresis loop. 
20 In still another embodiment, the comparing step is carried out by: (i) filtering the first 
and second QT- and RR- interval data sets by moving average smoothing; (ii) 
generating a smoothed hysteresis loop from the filtered first and second QT- and RR- 
interval data sets; and then (Hi) determining a measure of the domain inside the 
hysteresis loop. 

25 A further aspect of the present invention is a computer system for assessing 

cardiac ischemia in a subject to provide a measure of cardiac or cardiovascular health 

in that subject, the system comprising: 

(a) means for providing a first QT- and RR- interval data set from the subject 

during a stage of gradually increasing heart rate; 
30 (b) means for providing a second QT- and RR- interval data set from the 

subject during a stage of gradually decreasing heart rate; 
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■ 

(c) means for comparing the first QT- and RR- interval data set to the second 
QT- and RR- interval data set to determine the difference between the data sets; and 

i 

(d) means for generating from the comparison of step (c) a measure of cardiac 

* 

ischemia during stimulation in the subject, wherein a greater difference between the 

* 

5 first and second data sets indicates greater cardiac ischemia and lesser cardiac or 

■ 

cardiovascular health in the subject 

■ 

* 

A still further aspect of the present invention is a computer program product 

* * 

for assessing cardiac ischemia in a subject to provide a measure of cardiac or 
cardiovascular health in that subject, the computer program product comprising a 
10 . computer usable storage medium having computer readable program code means 
embodied in the medium, the computer readable program code means comprising: 

(a) computer readable program code means for comparing a first QT- and RR- 
interval data set to a second QT- and RR- interval data set to determine the difference 
between the data sets; and 

* 

15 (b) computer readable program code means for generating from the 

comparison of step (c) a measure of cardiac ischemia during stimulation in the 
subject, wherein a greater difference between the first and second data sets indicates 
greater cardiac ischemia and lesser cardiac or cardiovascular health in the subject. 

While the present invention is described herein primarily with reference to the 

20 use of QT- and RR- interval data sets, it will be appreciated that the invention may be 

. . implemented in simplified form with the use of RR- interval data sets alone. For use 
in the claims below, it will be understood that the term "RR-interval data set" is 
intended to be inclusive of both the embodiments of QT- and RR-interval data sets 
and RR-interval data sets alone, unless expressly subject to the proviso that the data 

25 set does not include a QT- interval data set. 

. The present invention is explained in greater detail in the drawings herein and 

■ 

the specification set forth below. 

Brief Description of the Drawings 
30 Figure 1 is a schematic graphic representation of the action potential in. 

cardiac muscle summed up over its volume and the induced electrocardiogram (ECG) * 
recorded on a human's body surface. 
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Figure 2A depicts the equations used in a simplified mathematical model of 
periodic excitation. 

Figure 2B depicts a periodic excitation wave (action potential, u, and 
instantaneous threshold, v, generated by computer using a simplified mathematical 
5 model, the equations of which are set forth in Figure 2A. 

Figure 2C depicts a family of four composite dispersion-restitution curves 
corresponding to four values of the medium excitation threshold. 

Figure 3 is a block diagram of an apparatus for carrying out the present 
method. 

10 Figure 4A is a block diagram of the .processing steps fctf data acquisition and 

analysis of the present invention. 

. Figure 4B is an alternative block diagram of the processing steps for data 
acquisition and analysis of the present invention. 

Figure 5 illustrates experimental QT-interval versus RR-interval hysteresis 
15 loops for two healthy male (23 year old, thick line and 47 year old, thin line) subjects 
plotted on the composite dispersion-restitution curve plane. 

Figure 6 provides examples of the QT-RR interval hysteresis for two male 

■ 

subjects, one with a conventional ECG ST-segment depression (thin line) and one 
with a history of a myocardial infarction 12 years prior to the test (thick line). The 

20 generation of the curves is explained in greater detail in the specification below. 

Figure 7 illustrates sensitivity of the present invention and shows two QT-RR 
interval hysteresis loops for a male subject, the first one (thick lines) corresponds to 
the initial test during which an ST-segment depression on a conventional ECG was 
observed, and the second one shown by thin lines measured after a period of regular 

25 exercise. 

« ■ 

Figure 8 illustrates a comparative cardiac ischemia analysis based on a 
particular example of a normalized measure of the hysteresis loop area. <CII> = (CII 
-Cllmini^CIImax - CH,™,) ("CII" means "cardiac ischemia index"). Oi, Xi, and Yi 
represent human subject data. Xj represents data collected from one subject (0.28 - 
30 0.35) in a series of tests (day/night testing, run/walk, about two months between 
tests); exercise peak heart rate ranged from 120 to 135. Yj represents data collected 
from one subject (0.46 - 0.86) in a series of tests (run/walk, six weeks between tests 
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* 

before and after a period of regular exercise stage); exercise peak heart rate ranged 
from 122 to 146. Black bars indicate a zone (<CH> less than 0.70) in which a 
conventional ST depression method does not detect cardiac ischemia. The 
conventional method may detect cardiac ischemia only in a significantly narrower 
5 range indicated by high white bars (Y 2 , Y3, O7: <CII> greater than 0.70). 

Figure 9 illustrates a typical rapid peripheral nervous system and hormonal 
control adjustment of the QT and RR interval to an abrupt stop in exercise (that is, an 
abrupt initiation of a rest stage). 

Figure 10 illustrates a typical slow (quasi-stationary) QT and RR interval 
10 adjustment measured during gradually increasing and gradually decreasing cardiac 
stimulation. 

Figure 11 demonstrates a block-diagram of the data processing by the method 
of optimized consolidation of a moving average, exponential and polynomial fitting 
(Example 10, steps 1-8). 
15 Figure 12 demonstrates results of the processing throughout steps 1 to 8 of 

Example 10. Upper panels show QT and RR data sets processed from steps 1 to 3 . 

* * 

• 4 f 

. (from left to right respectively), and the QT/RR hysteresis loop after step 1. The 
exponential fitting curves (step 3) are shown in gray in the first two panels. Low panels 

* 

show the same smooth dependencies after processing from step 4 to a final step 8. Here 
20 the CII (see right . low . . panel) is equal to a ratio 

S/{((T^(t end -t»)-7Z(WW„ -^-r* r (0))}(see examplelO, section 7). 

Figure 13 demonstrates a block-diagram of the data processing by the method 
of a sequential moving average (Example 11, steps 1-3) 

Figure 14 demonstrates results of the processing throughout steps 1 to 2 of 
25 Example 11. Upper panels show processed QT and RR data sets, and the QT/RR 
hysteresis loop after step 1 (from left to right respectively). Low panels show the same 
smooth dependencies after the second moving average processing and a final step 3. 
Figure 15 shows a general data flow chart for major steps in optimized 
, nonlinear transformation method. The left-hand side aind the right-hand side boxes 
30 describe similar processing stages for the RR and QT-intervals, respectively. 

Figure 16 shows a detailed data flow chart for one data subset, {*\Z*rr} or 

* m 

{A^qt} during stages. (1 a/b) through (3 a/b) in Figure 15. The preliminary stage, 
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boxes 1 through 7, uses a combination of traditional data processing methods and 
includes: moving averaging (1), determination of a minimum region (2), fitting a 
quadratic parabola to the data in this region (3), checking consistency of the result (4), 
finding the minimum and centering data at the minimum (5) and (6), conditionally 

5 sorting the data (7). Stages (8) through (11) are based on the dual-nonlinear' 

■ 

transformation method for the non-linear regression. 

i 

Figure 17 displays the nonlinear transformation of a filtered RR-interval data 
set. Panel A shows the data on the original (/,y)-plane, the minimum is marked with 

■ 

an asterisk inside a circle. Panels B and C show the transformed sets on the (f,u)- 

- § 

1 0 plane, for j = 1 and for j = 2, respectively. The image of the minimum -is also marked 
with an encircled asterisk. Note that the transformed data sets concentrate around a 
monotonously growing (average) curve with a clearly linear portion in the middle. 
Figure 18 is similar to Figure 17 but for a QT-interval data set. 
Figure 19 displays an appropriately scaled representation for the family of 
15 functions i;(a,P,T) for fifteen values of parameter P varying with the step Ap^O.l from 
P=-0.9 (lower curve) through p=0 (medium, bold curve), to P=0.5 (upper curve). The 
function ^(a,P,x) is continuous in all three variables and as a function of t has a unit 
slope at x=0, £(a,P,0)=l, 

Figure 20 shows an example of full processing of the RR and QT data sets for 
20 one patient. Panels A and C represent RR and QT data sets and their fit. Panel B shows 
the corresponding ascending and descending curves and closing line on the (Trr,Tqt)- 
plane, on which the area of such a hysteresis loop has the dimension of time-squared. 
Panel D shows a hysteresis loop on the (/^TQ^-plane, where Jrr=1/Trr is the heart 
rate, on which the loop area is dimensionless. The total error-for Panel A is 2.2% and 
25 for Panel C is 0.8%. 

Detailed Description of the Preferred Embodiments 

The present invention is explained in greater detail below. This description is 
not intended to be a detailed catalog of all the different manners in which particular 
30 elements of the invention can be- implemented, and numerous variations will be 
apparent to those skilled in the art based upon the instant disclosure. 
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As will be appreciated by one of skill in the art, certain aspects of the present 
invention may be embodied as a method, data processing system, or computer 
program product. Accordingly, certain aspects of the present invention may take the 
form of an entirely hardware embodiment, an entirely software embodiment, or an 
5 embodiment combining software and hardware aspects. Furthermore, certain aspects 
of the present invention may take the form of a computer program product on a 

■ * 

■ 

computer-usable storage medium having computer readable program code means 
embodied in the medium. Any suitable computer readable medium may be utilized 

» ■ 

including, but not limited to, hard disks, CD-ROMs, optical storage devices, and 

10 magnetic storage devices. 

Certain aspects of the present invention are described below with reference to 
flowchart illustrations of methods, apparatus (systems), and computer program 
products. It will be understood that each block of the flowchart illustrations, and 
combinations of blocks in the flowchart illustrations, can be implemented by 

15 computer program instructions. These computer program instructions may be 
provided to a processor of a general purpose computer, special purpose computer, or 
other programmable data processing apparatus to produce a machine, such that the 
instructions, which execute via the processor of the computer of other programmable 
data processing apparatus, create means for implementing the functions specified in 

i 

20 the flowchart block or blocks. 

Computer program instructions, may also be stored in a computer-readable 

• * . 

■ i 

* 4 # 

memory that can direct a computer or other programmable data processing apparatus 
to function in a particular manner, such that the instructions stored in the computer- 
readable memory produce an article of manufacture including instruction means 
. 25 which implement the function specified in the flowchart blpck or blocks. 

Computer program instructions may also be loaded onto a computer or other 
programmable data processing apparatus to cause a series of operational steps to be 
performed on the computer or other programmable apparatus to produce a computer 
implemented process such that the instructions which execute on the computer or 
30 other programmable apparatus provide steps for implementing the functions specified , 
in the flowchart block or blocks. 
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1. Definitions. 

"Cardiac ischemia" refers to a lack of or insufficient blood supply to an area of 
cardiac muscle. Cardiac ischemia usually occurs in the presence of arteriosclerotic 
occlusion of a single or a group of coronary arteries. Arteriosclerosis is a product of a 
5 lipid deposition process resulting in fibro-fatty accumulations, or plaques, which grow 

■ 

* 

on the internal walls of coronary arteries. Such an occlusion compromises blood flow 
through the artery, which reduction then impairs oxygen supply to the surrounding 
tissues during increased physiological need — for instance, during increased exercise 
loads. In the later stages of cardiac ischemia (e.g., significant coronary artery 
10 occlusion), the blood supply may be insufficient even while the cardiac muscle is at 
rest. However, in its earlier stages such ischemia is reversible in a manner analogous 
to how the cardiac muscle is restored to normal function when the oxygen supply to it 

■ 

returns to a normal physiological level. Thus, ischemia that may be detected by the 
present invention includes episodic, chronic and acute ischemia. 
15 "Exercise" as used herein refers to voluntary skeletal muscle activity of a 

subject that increases heart rate above that found at a sustained stationary resting state. 
Examples of exercise include, but are not limited to, cycling, rowing, weight-lifting, 

-» 

walking, running, stair-stepping, etc., which may be implemented on a stationary 
device such as a treadmill or in a non-stationary environment. 
20 "Exercise load" or "load level" refers to the relative strenuousness of a 

particular exercise, with greater loads or load levels for a given exercise producing a 
greater heart rate in a subject. For example, load may be increased in weight-lifting 
by increasing the amount of weight; load may be increased in walking or running by 
increasing the speed and/or increasing the slope or incline of the walking or running 

» 

25 surface; etc. 

■ 

"Gradually increasing" and "gradually decreasing" an exercise load refers to 
exercise in which the subject is caused to perform an exercise under a plurality of 
different sequentially increasing or sequentially decreasing loads. The number of 
steps in the sequence can be infinite so the terms gradually increasing and gradually 
30 decreasing loads include continuous load increase and decrease, respectively. 

"Intervening rest", when used to refer to a stage follpwing increased cardiac 
stimulation, refers to a stage of time initiated by a sufficiently abrupt decrease in heart 
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stimulation {e.g., an abrupt decrease in exercise load) so that it evokes a clear 
sympatho-adrenal response. Thus, an intervening rest stage is characterized by a 

r 

rapid sympatho-adrenal adjustment (as further described in Example 8 below), and the 
. inclusion of an intervening rest stage precludes the use of a quasi-stationary exercise 
5 (or stimulation) protocol (as further described in Example 9 below). 

-"Hysteresis" refers to a lagging of the physiological effect when the external 
conditions are changed. ... 

"Hysteresis curves" refer to a pair of curves in which one curve reflects the 
response of a system to a first sequence of conditions, such as gradually increasing 
10 heart rate, and the other curve reflects the response of a system to a second sequence 

» 

of conditions, such as gradually decreasing heart rate. Here both sets of conditions are 

* 

essentially the same--/.*., consist of the same (or approximately the same) steps— but 
are passed in different order in the course of time. A "hysteresis loop" refers to a loop 
formed by the two contiguous curves of the pair. 
15 "Electrocardiogram" or "ECG" refers to a continuous or sequential record (or 

a set of such records) of a local electrical potential field obtained from one or more 
locations outside the cardiac muscle. This field is generated by the combined 

* > 

electrical activity (action potential generation) of multiple cardiac' cells. The 
recording electrodes may be either subcutaneously implanted or may be temporarily 
20 attached to the surface of the skin of the subject, usually in the thoracic region. An 
ECG record typically includes the single-lead ECG signal that represents a potential 
difference between any two of the recording sites including the site with a zero or 
ground potential. 

"Quasi-stationary conditions" refer to a gradual change in the external 
25 conditions and/or the physiological response it causes that occurs much slower than 
any corresponding adjustment due to sympathetic/parasympathetic and hormonal 
control. If the representative time of the external conditions variation is denoted by 
Text, and Xi n t is a representative time of the fastest of the internal, 
sympathetic/parasympathetic and hormonal control, then "quasi-stationary 
30 conditions" indicates r ext » ti nt {e.g., T ext is at least about five times greater than Tj nt ). 
"An abrupt change" refers to an opposite situation corresponding to a sufficiently fast 
change in the external conditions as compared with the rate 
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sympathetic/parasympathetic and hormonal control — that is, it requires that x ex t « Tint 
(e.g., Text is at least about five times less than T int ). In particular, "an abrupt stop" refers 
to a fast removal of the exercise load that occurs during time shorter than Tj nt ~ 20 or 
. 30 seconds (see Figure 9 below and comments therein). 
5 . "QT- and RR- data set" refers to a record of the time course of an electrical 

signal comprising action potentials spreading through cardiac muscle. Any single lead 
ECG record incorporates a group of three consecutive sharp deflections usually called - 
a QRS complex and generated by the propagation of the action potential's front 
through the ventricles. In contrast, the electrical recovery of ventricular tissue is seen . 

■ 

10 on the ECG as a relatively small deflection known as the T wave. The time interval 
between the cardiac cycles (i.e., between the maxima of the consecutive R-waves) is 
called an RJR-interval, while the action potential duration (i.e., the time between the 
beginning of a QRS complex and the end of the ensuing T-wave) is called a QT- 
interval. Alternative definitions of these intervals can be equivalently used in the 

1 5 framework of the present invention. For example, an RR-interval can be defined as 

* 

the time between any two similar points, such as the similar inflection points, on two 
consecutive R-waves, or any other way to measure cardiac cycle length. A QT- 
interval can be defined as the time interval between the peak of the Q-wave and the 
peak of the T wave. It can also be defined as the time interval between the beginning 
20 (or the center) of the Q-wave and the end of the ensuing T-wave defined as the point 
on the time axis (the base line) at which it intersects with the linear extrapolation of 
the T-wave's falling branch and started from its inflection point, or any other way to 

► 

measure action potential duration. An ordered set of such interval durations 
simultaneously with the time instants of their beginnings or ends which are 

■ 

25 accumulated on a beat to beat basis or on any given beat sampling rate basis form a 
corresponding QT- and RR-interval data set. Thus, a QT- and RR- interval data set 
will contain two QT-interval related sequences {TQT,i,TQT,2v">TQT, n >} and {ti,t2,...,tn}, 
and will also contain two RR-interval related sequences {Trr^Trr^.^Trr^J and 
{ti,t2,...,t n } (the sequence {ti,t2,...,t n } may or may not exactly coincide with the similar 

30 sequence in the QT data set). 

In the following definitions, C[a,b] shall denote a set of continuous functions 
fit) on a segment [a,b]. {t\}, i=l,2,..., N 9 denotes a set of points from [a,b], i.e. 
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{/j}={*i: a<t\ < b f /=1,2,...^V} and {/ft)}> where feC[a,b], denotes a set of values of 
the function / at the points {t\}. In matrix operations the quantities v={t\}, ^{/fa)}, 
are treated as column vectors. £n shall denote a iV-dimensional metric space with the 

metric Rn(x,y), x,yGEn. (Ru(x,y) is said to be a distance between points x and y.) A 

* 

b 

(total) variation \f [F\ is defined for any absolutely continuous function F from 

a * 

* 

C[a,b] as the integral (a Stieltjes integral) 

b b b 



V [F(/)]^ \\dF{t)\=\\F{t)\dt. (D.l) 



For a function F monotonic on segment [a,b] its variation is simply |F(a)-F(b)|. If a 
function F(t) has alternating maxima and minima, then the total variation of F is the 
1 0 sum of its variations on the intervals of monotonicity. For example, if the points of 

■ 

minima and maxima are x\=a, * 2 , x\c=b then 

V[F{t)] = f,\F(x,)-F(x M )\. (D.2) 

a M 



Fitting (best fitting): Let C[a,b] be a subset of C[a,b]. A continuous 
15 function f[t)>fe C[a,b] is called the (best) fit (or the best fitting) function of class 
C [a,b] with respect to metric i? N to a data set {x-Ji} (i=l ,2, . . . , N} if 



M{Kt\)},{Xi})= min (D.3) 
The minimum value of i?N is then called the error of the fit. The functions/^) from 

• ■ 

20 C [a,b] will be called *na/ functions. 

In most cases £n is implied to be an Euclidean space with an Euclidean metric. 
The error R^ then becomes the familiar mean-root-square error. The fit is performed 

on a subset. C[a,b] since it usually implies a specific parametrization of the trial 
functions and/or such constrains as the requirements that the -trial functions pass 

* * 

25 through a given point and/or have a given value of the slope at a given point. 
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A smoother function (comparison of smoothness): Let fit) and g(f) be 
functions from C[a,b] that have absolutely continuous derivatives on this segment. 
The function J{t) is smoother than the function g(t) if 



VW0]<Vte(O], . (D.4) 



a 



and 



W(01<W(OL (D.5) 



a 



. where the prime denotes a time derivative, and a strict inequality holds in at least one 
of relations (D.4) and (D.5). 
10 A smoother set: A set {xJi} (j-1,2,..., N} is smoother than the set {x),t)} 

(/-1,2,..., N} if the former can be fit with a smoother function /(/) of the same class 

within the same or smaller error than the latter. 

» 

Smoothing of a data set: A (linear) transformation of a data set (x,t)={x\,t\) 
0—1,2,..., No} into another set (y/r)={yj,Tj} (/=1,2,..., N\} of the form 



y=A-x, t=B% (D.6) 



where A and B are N\ x No matrices, is called a smoothing if the latter set is smoother 
than the former. One can refer to {yj,Tj} as a smoothed set 
20 A measure of a closed domain: Let £1 be a singly connected domain on the 

plane (t,1) with the boundary formed by a simple (i.e., without self-intersections) 
continuous curve. A measure M of such a domain Q. on the plane (t, T) is defined as 
the Riemann integral 

; M = \\p{r,T)dxdT . (D.7) 
n 

25 where p(x,7) is a nonnegative (weight) function on £1 

Note that when p(t,7)=l the measure M of the domain coincides with its area, 
A; when p(t,7)=1/t 2 , the measure, M, has the meaning of the area, A\ of the domain 
Q' on .the transformed plane (f,T), where /=7A can be understood as the heart rate 
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. since the quantity x has the meaning of RR-interval. [The domain CI' is the image of 
domain CI under the mapping (x,7)— ►(1/t,7).] 

i 

2. Dispersion/restitution curves. 
5 Figure 1 illustrates the correspondence between the temporal phases of the 

periodic action potential (AP, upper graph, 20) generated inside cardiac muscle and 
summed up over its entire volume and the electrical signal produced on the body 
surface and recorded as an electrocardiogram (ECG, lower graph, 21). The figure 
depicts two regular cardiac cycles. During the upstroke of the action potential the 
10 . QRS-complex is formed. It consists of three waves, Q 7 R, and S, which are marked on 
the lower panel. The recovery stage of the iaction potential is characterized by its fall 
off on the AP plot and by the T-wave on the ECG plot. One can see that the action 
potential duration is well represented by the time between Q and T waves and is 

■ 

conventionally defined as the QT interval, measured from the beginning of the Q 
15 wave to the end of the following T wave. The time between consecutive R-waves (RR 
interval) represents the duration of a cardiac cycle, while its reciprocal value 
represents the corresponding instantaneous heart rate. 

Figure 2 illustrates .major aspects of the process of propagation of a periodic 
action potential through cardiac tissue and the formation of a corresponding 
20 composite dispersion-restitution curve. The tissue can be considered as a continuous 
medium and the propagation process as a repetition at each, medium point of the 
. consecutive phases of excitation and recovery. The former phase is characterized by a 

fast growth of the local membrane potential (depolarization), and the latter by its 

« 

^ 

return to a negative resting value (repolarization). The excitation phase involves a 

* 

25 very fast (-O.lrns) decrease in the excitation threshold and the following development 

« 

of a fast inward sodium current that causes the upstroke of the action potential 

■ ■ . » 

(~lms). Next, during an intermediate plateau phase (~200ms) sodium current is 
inactivated; calcium and potassium currents are developing while the membrane is 
temporarily unexcitable (i.e., the threshold is high). During the next recovery phase 

30 (~100ms), a potassium current repolarizes the membrane so it again becomes 

i 

excitable (the excitation threshold is lowered). 
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The complicated description of a multitude of ionic currents involved in the 
process can be circumvented if one treats the process directly in terms of the local 

■ 

membrane potential, u, and a local excitation threshold, v. Such a mathematical 
description referred to as the CSC model, was developed by Chernyak, Stairobin, & 
5 Cohen (Phys. Rev. Lett., 80, pp;5675-5678, 1998) and is presented as a set of two 
Reaction-Diffusion (RD) equations in panel A. The left-hand side of the first equation 
describes local accumulation of the electric charge on the membrane, the first term in 
the right-hand side describes Ohmic coupling between neighboring points of the 
medium, and the term i(u,v) represents the transmembrane current as a function of the 

10 membrane potential and the varying excitation threshold (e is a small constant, the 
ratio of the slow recovery rate to the fast excitation rate). A periodic solution (a wave- 
train) can be found analytically for some particular functions i(u,v) and g{u,v). The 
wave-train shown in panel B has been calculated for g(w,v)=^w-H> r -v, where C, and v r 
are appropriately chosen constants (v r has the meaning of the initial excitation 

15 threshold and is the main determinant of the medium excitability). The function i(u f v) 
was chosen to consist of two linear pieces, one for the sub-threshold region, w<v, and 
one for supra-threshold region, u>v. That is i(u,v)=X& when u<v, and i(w,v)=k e k(«- 
i/ex) when w>v, where X T and Xex are membrane chord conductances in the resting 
(tt=0) and excited (u=u cx ) states, respectively. The resting state u-0 is taken as the 

■ 

20 origin of the potential scale. We used such units that A* x =l and u ex =l. (For details see 
Chernyak & Starobin, Critical Reviews in Biomed. Eng. 27, 359-414 (1999)). 

A medium with higher excitability, corresponding to the tissue with better 

* 

conduction, gives rise to a faster, more robust action potential with a longer Action 
Potential Duration (APD). This condition also means that a longer-lasting excitation 
25 propagates- faster. Similarly, a wave train with a higher frequency propagates slower 

4 

since the medium has less time to recover from the preceding excitation and thus has 
a lower effective excitability. These are quite generic features that are incorporated in 
the CSC model. In physics, the relation between the wave's speed, c, and its 
frequency,/, or its period, T=J/f, is called a dispersion relation. In the CSC model the 
30 - dispersion relation can be obtained in an explicit form" T-Fj{c), where Fj is a known 
function of c and the medium parameters. The CSC model also allows us to find a 
relation between the propagation speed and the APD, Tap, in the explicit form 
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^ap-^ap(c), which represents the restitution properties of the medium. In the medical 

* . 

literature, the restitution curve is Tap versus diastolic interval Jbi, which differently • 
makes a quite similar physical statement. One can consider a pair of dispersion and 
restitution relations {T=Ft(c),Tap=Fap(c)} as a parametric representation of a single 

* < * • 

5 curve on the (T, 7Ap)-plane as shown in panel C (Figure 2). Such a curve (relation) 
shall be referred to as a composite dispersion-restitution curve (relation) and can be 

■ * 

directly obtained from an experimental ECG recording by determining the QT-'RR 
interval data set and plotting 7q T versus Jrr. A condition that the experimental 
{JqtjTrr} data set indeed represents the composite dispersion-restitution relation is 
10 the requirement that the data ' are collected under quasi-stationary conditions. 
Understanding this fact is a key discovery for the present invention. - 

* 

* 

3. Testing methods. 

The methods of the present invention are primarily intended for the testing of 
15 human subjects. Virtually any human subject can be tested by the methods of the 
present invention, including male, female, juvenile, adolescent, adult, and geriatric 
subjects. The methods may be carried out as an initial screening test on subjects for 
which no substantial previous history or record is available, or may be carried out on a 
repeated basis on the same subject (particularly where a comparative quantitative 

• - 

20 indicium of an individual's cardiac health over time is desired) to assess the effect or 
influence of intervening events and/or intervening therapy on that subject between 
testing sessions. 

■ 

As noted above, the method of the present invention generally comprises (a) 
collecting a first QT- and RR- interval data set from said subject during a stage of 

■ 

25 gradually increasing heart rate; (b) collecting a second QT- and RR- interval data set 
from said subject during a stage of gradually decreasing heart rate; (c) comparing said 

■ 

first QT- and RR- interval data set to said second QT- and RR- interval data set to 
determine the difference between said data sets; and (d) generating from said 
comparison of step (c) a measure of cardiac ischemia in the subject. A greater 
30 difference between the first and second data sets indicates greater cardiac ischemia 
and lesser cardiac or cardiovascular health in that subject. 



WO 02/00113 PCTYUS01/20306 

» 

■ 

-23- 

The stages of gradually increasing and gradually decreasing heart rate are 
carried out in a manner that maintains during both periods essentially or substantially 
the same stimulation of the heart by the peripheral . nervous and hormonal control 
systems, so that it is the effect of cardiac ischemia rather than that of the external 
5 control which is measured by means of the present invention. This methodology can 
be carried out by a variety of techniques, with the technique of conducting two 
consecutive stages of gradually increasing and gradually decreasing exercise loads (or 
average heart rates) being currently preferred. 

The stage of gradually increasing exercise load (or increased average heart 
10 rate) and the stage of gradually decreasing exercise load (or decreased average heart 
rate) may be the same in duration or may be different in duration. In general, each 
. stage is at least 3, 5, 8, or 10 minutes or more in duration. Together, the duration of 
the two stages may be from about 6, 10, 16 or 20 minutes in duration to about 30, 40, 
or 60 minutes in duration or more. The two stages are preferably carried out 
15 sequentially in time — that is, with one stage following after the other substantially 
immediately, without an intervening rest stage. In the alternative, the two stages may 
be carried out separately in time, with an intervening "plateau" stage {e.g., of from 1 
to 5 minutes) during which cardiac stimulation or exercise load is held substantially 

■ 

constant, before the -stage of decreasing load is initiated. 
20 The exercise protocol may include the same or different sets of load steps 

during the stages of increasing or decreasing heart rates. For example, the peak load in 
each stage may be the same or different, and the minimum load in each stage may be 
the same or different. In general, each stage consists .of at least two or three different 
load levels, in ascending or descending order depending upon the stage. Relatively 

4 

25 high load levels; which result in relatively high heart rates, can be used but are not 
essential. An advantage of the present invention is that its sensitivity allows both 
exercise procedures to be carried out at relatively low load levels that do not unduly 

* 

increase the pulse rate of the subject. For example, the method may be carried out so 
that the heart rate of the subject during either the ascending or descending stage (or 
30 both) does not exceed about 140, 120, or even 100 beats per minute, depending upon 
the condition of the subject. Of course, data collected at heart rates above 100, 120, 
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or 140 beats per minute may also be utilized if desired, again depending upon the 
condition of the subject. 

For example, for an athletic or trained subject, for the first or ascending stage, 
a first load level may be selected to require a power output of 60 to 100 or 150 watts 
5 by the subject; an intermediate load level may be selected to require a power output of 
100 to 150 or 200 watts by the subject; and a third load level may be selected to 
require a power output of 200 to 300 or 450 watts or more by the subject. For the 
second or descending stage, a first load level may be selected to require a power 
output of 200 to 300 or 450 watts or more by the subject; an intermediate or second 

10 load level may be selected to require a power output of 100 to 1 50 or 200 watts by the 
subject; and a third load level may be selected to require a power output of 60 to 100 
or 150 watts by the subject. Additional load levels may be included before, after, or 
between all of the foregoing load levels as desired, and adjustment between load 
levels can be carried out in any suitable manner, including step- wise or continuously. 

15 In a further example, for an average subject or a subject with a history of 

cardiovascular disease, for the first or ascending stage, a first load level may be 
selected to require a power output of 40 to 75 or 100 watts by the subject; an 
intermediate load level may be selected to require a power output of 75 to 100 or 150 
watts by the subject; and a third load level may be selected to require a power output 

20. of 125 to 200 or 300 watts or more by the subject. For the second or descending 
stage, a first load level may be. selected to require a power output of 125 to 200 or 300 
watts or more by the subject; an intermediate or second Toad level may be selected to 
require a power output of 75 to 100 or 150 watts by the.subject; and a third load level 

< 

may be selected to require a power output of 40 to 75 or 100 watts by the subject. As 
25 before, additional load levels may be included before, after, or between all of the 

foregoing load levels as desired, and adjustment between load levels can be carried 

. <■ 

out in any suitable manner, including step-wise or continuously. 

The heart rate may be gradually increased and gradually decreased by 
subjecting the patient to a predetermined schedule of stimulation. For example, the 
30 patient may be subjected to a gradually increasing exercise load and gradually 
decreasing exercise load, or gradually increasing electrical or pharmacological 
stimulation and gradually decreasing electrical or pharmacological stimulation, 



I* 



WO 02/00113 



PCT/USO 1/20306 



-25- 

according to a predetermined program or schedule. Such a predetermined schedule is 

* 

without feedback of actual heart rate from the patient. In the alternative, the heart rate 
of the patient may be gradually increased and gradually decreased in response to 
actual heart rate data collected from concurrent monitoring of said patient. Such a 
5 system is a feedback system. For example, the heart rate of the patient may be 
monitored during the test and the exercise load (speed and/or incline, in the case of a 

> 

treadmill) can be adjusted so that the heart rate varies in a prescribed way during both 
stages of the test. The monitoring and control of the load can be accomplished by a 
computer or other control system using a simple control program and an output panel 

10 connected to the control system and to the exercise device that generates an analog 
signal to the exercise device. One advantage of such a feedback system is that* (if 
desired) the control system can insure that the heart rate increases substantially 
linearly during the first stage and decreases substantially linearly during the second 
stage. ' 

1 5 The generating step (d) may be carried out by any suitable means, such as by 

generating curves from the data sets (with or without actually displaying the curves), 
and then (i) directly or indirectly evaluating a measure (e.g., as defined in the integral 
theory) of the domain (e.g., area) between the hysteresis curves, a greater measure 
indicating greater cardiac ischemia in said subject, (ii) directly or indirectly 

20 comparing the shapes (e.g., slopes or derivatives thereof) of the curves, with a greater 
difference in shape indicating greater cardiac ischemia in the subject; or (Hi) 
combinations of (i) and (ii). Specific examples are given in Example 4 below. 

The method of the invention may further comprise the steps of (e) comparing 
the measure of cardiac ischemia during exercise to at least one reference value (e.g., a 

25 mean, median or mode for the quantitative indicia from a population or subpopulation 
of individuals) and then (f) generating from the comparison of step (e) at least one 
quantitative indicium of cardiovascular health for said subject. Any such quantitative 
indicium may be generated on a one-time basis (e.g., for assessing the likelihood that 
the subject is at risk to experience a future ischemia-related cardiac incident such as 

30 myocardial infarction or ventricular tachycardia), or may be generated to monitor the 
progress of the subject over time, either in response to a particular prescribed 

■ 

cardiovascular therapy, or simply as an ongoing monitoring of the physical condition 
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of the subject for improvement or decline (again, specific examples are given in 

r 

Example 4 below). In such a case, steps (a) through (j) above are repeated on at least 
one separate occasion to assess the efficacy of the cardiovascular therapy or the 
progress of the subject. A decrease in the difference between said data sets from 

* ■ 
- 

5 . before said therapy to after said therapy, or over time, indicates an improvement in 
cardiac health in said subject from said cardiovascular therapy. Any suitable 
cardiovascular therapy can be administered, including but not limited to, aerobic 
exercise, muscle strength building, change in diet, nutritional supplement, weight loss, 
smoking cessation, stress reduction, pharmaceutical treatment (including gene 

10 therapy), surgical treatment (including both open heart and closed heart procedures 
such as bypass, balloon angioplasty, catheter ablation, etc.) and combinations thereof 
The therapy or therapeutic intervention may be one that is approved or one 
that is experimental. In the latter case, the present invention may be implemented in 
the context of a clinical trial of the experimental therapy, with testing being carried 

15 out before and after therapy (and/or during therapy) as an aid in determining the 

* • % 

efficacy of the proposed therapy. 

4. Testing apparatus. 

Figure 3 provides an example of the apparatus for data acquisition, processing 

20 and analysis by the present invention. Electrocardiograms are recorded by an ECG 
recorder 30, via electrical leads placed on a subject's body. The ECG recorder may 
be, for example, a standard multi-lead Holter recorder or any other appropriate 
recorder. The analog/digital converter 31 digitizes the signals recorded by the ECG 
recorder and transfers them to a personal computer 32, or other computer or central 

25 processing unit, through a standard external input/output port. The digitized ECG data 
can then be processed by standard computer-based waveform analyzer software. 
Composite dispersion-restitution curves and a cardiac or cardiovascular health 
indicium or other quantitative measure of .the presence, absence or degree of cardiac 
ischemia can then be calculated automatically in the computer through a program 

30 (e.g., Basic, Fortran, C++, eta) implemented therein as software, hardware, or both 

> ■ 

hardware and software. 



WO 02/00113 



PCT/US01/20306 



-27- 

■ 

Figure 4A and Figure 4B illustrate the major steps of digitized data 
processing in order to generate an analysis of a QT-RR data set collected from a * 
subject during there-and-back quasi-stationary changes in physiological conditions. 
The first four steps in Figure 4 A and Figure 4B are substantially the same. The 
5 digitized data collected from a multi-lead recorder are stored in a computer memory 

* • 

for each lead as a data array 40a, 40b. The size of each data array is determined by the 
durations of the ascending and descending heart rate stages and a sampling rate used 
by the waveform analyzer, which processes an incoming digitized ECG signal. The 
waveform analyzer software first detects major characteristic waves (Q,R,S and T 
10 waves) of the ECG signal in each particular lead 41a, 41b. Then in each ECG lead it 
determines the time intervals between consecutive R waves and the beginning of Q 

i 

r ■ 

and the end of T waves 42a, 42b. Using these reference points it calculates heart rate 
and RR- and QT- intervals. Then, the application part of the software sorts the 
intervals for the ascending and descending heart rate stages 43a, 43b. The next two 

* 

15 steps can be made in one of the two alternative ways shown in Figures 4 A and 4B, 

k 

respectively. The fifth step as shown in Figure 4A consists of displaying by the 
application part of software QT- intervals versus RR- intervals 44a, separately for the 

■ 

ascending and descending heart rate stages effected by there-and-back gradual 
changes in physiological conditions such as exercise, pharmacological/electrical 

20 stimulation, etc.. The same part of the software performs the next step 45a, which is 
smoothing, filtering or data fitting, using exponential or any other suitable functions, 
in order to obtain a sufficiently smooth curve Tqt=F(Trr) for each stage. An 
alternative for the last two steps shown in Figure 4B requires that the application part 
of the software first averages, and/or filters and/or fits, using exponential or any other 

25 suitable functions, the QT intervals as functions of time for both stages and similarly 
processes the RR-interval data set to produce two sufficiently smooth curves 
Jqt^qtCO m & ^rr = ^rr(0> ea °h including the ascending and descending heart rate 
branches 44b. At the next step 45b the application part of the software uses this 
parametric representation to eliminate time and generate and plot a sufficiently 

30 smooth hysteresis loop Tqt=F(Trr)' The following steps shown in Figures 4A ' and 

* 

Figure 4B are again substantially the same. The next step 46a, 46b performed by the 
application part of the software can be graphically presented as closing the two branch 
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hysteresis loop with an appropriate interconnecting or partially connecting line, such 
as a vertical straight line or a line connecting the initial and final points, in order to 
produce a closed hysteresis loop on the (7qt, JRR)-plane. At the next step 47a, 47b the 
application software evaluates for each ECG lead an appropriate measure of the 
5 domain inside the closed hysteresis loop. A measure, as defined in mathematical 
integral theory, is a generalization of the concept of an area and may include 
appropriate weight functions increasing or decreasing the contribution of different 
portions of the domain into said measure. The final step 48a, 48b of the data 
processing for each ECG lead is that the application software calculates indexes by 

10 appropriately renormalizing the said measure or any monotonous functions of said 
measure. The measure itself along with the indexes may reflect both the severity of 
the exercise-induced ischemia, as well as a predisposition to local ischemia that can 
be reflected in some particularities of the shape of the measured composite dispersion- 
restitution curves. The results of all above-mentioned signal processing steps may be 

15 used to quantitatively assess cardiac ischemia and, as a simultaneous option, 
cardiovascular system health of a particular individual under test. 

Instead of using the (Tqt, 7RR)-plane a similar data processing procedure can 
equivalently be performed on any plane obtained by a non-degenerate transformation 
of the (rQ T ,7RR)-plane, such as (7q T ,./rk) where fysrl/T^ is the heart rate or the like. 

20 Such a transformation can be partly or fully incorporated in the appropriate definition 

« 

• of the said measure. 

The present invention is explained in greater detail in . the non-limiting 
examples set forth below. ' 



25 EXAMPLE 1 

( * 

Testing Apparatus 

A testing apparatus consistent with Figure 3 was assembled. The 
electrocardiograms : are recorded by an RZ152PM12 Digital ECG Holter Recorder 
(ROZINN ELECTRONICS, INC.; 71-22 Myrtle Av, Glendale, New York, USA 
30 11385-7254), via 12 electrical leads with Lead-Lok Holter/Stress Test Electrodes 

i 

LL510 (LEAD-LOK, INC.; 500 Airport Way, P.CBox L, Sandpoint, ID, USA 
83864) placed on a subject's body in accordance with the manufacturer's instructions. 
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Digital ECG data are transferred .to a personal computer (Dell Dimension XPS 
TSOOMHz/Windows 98) using a 40 MB flash card (RZFC40) with a PC 700 flash 
card reader, both from Rozinn Electronics, Inc. Holter for Windows (4.0.25) 
waveform analysis software is installed in the computer, which is used to process data 
5 by a standard computer based waveform analyzer software. Composite dispersion- 
restitution curves and an indicium that provides a quantitative characteristic of the 
extent of cardiac ischemia are then calculated manually or automatically in the 
computer through a program implemented in Fortran 90. 

Experimental data were collected during an exercise protocol programmed in a 

1 0 Landice L7 Executive Treadmill (Landice Treadmills; 1 1 1 Canfield Av., Randolph, NJ 
07869). The programmed protocol included 20 step-wise intervals of a constant 
exercise load from 48 seconds to 1.5 minutes each in duration. Altogether these 
intervals formed two equal-in-duration gradually increasing and gradually decreasing 
exercise load stages, with total duration varying from 16 to 30 minutes. For each stage 

15 a treadmill belt speed and elevation varied there-and-back, depending on the subject's 
age and health conditions, from 1 .5 miles per hour to 5.5 miles per hour and from one 
to ten degrees of treadmill elevation, respectively. . 

EXAMPLES 2-6 
20 Human Hysteresis Curve Studies 

► * 

These examples illustrate quasi-stationary ischemia-induced QT-RR interval 
hystereses in a variety of different human subjects. These data demonstrate a high 

* 

sensitivity and the high resolution of the method. 

« 

25 EXAMPLES 2-3 

Hysteresis Curves in Healthy 
Male Subjects of Different Ages 
These examples were carried out on two male subjects with an apparatus and 
procedure as described in Example 1 above. Refemng to Figure 5, one can readily 
30 . see a significant difference in areas of hystereses between two generally healthy male 
subjects of different ages. These subjects (23 and 47 years old) exercised on a 

> 

treadmill according to a quasi-stationary 30-minute protocol with gradually increasing 
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and gradually decreasing exercise load. Here squares and circles (thick line) indicate a 
hysteresis loop for the 23 year old subject, and diamonds and triangles (thin line) 
correspond to a larger loop for the 47 year old subject. Fitting curves are obtained 
using the third-order polynomial functions. A beat sampling rate with which a 
5 waveform analyzer determines QT and RR intervals is equal to one sample per 
minute. Neither of the subjects had a conventional ischemia-induced depression of the 

» 

ECG-ST segments. However, the method of the present invention allows one to 
observe ischemia-induced hystereses that provide a satisfactory resolution within a 
conventionally sub-threshold range of ischemic events and allows one to 
10 quantitatively differentiate between the hystereses of the two subjects. 

EXAMPLES 4-5 

■ 

■ 

Hysteresis Curves for Subjects with ST 
Segment Depression or Prior Cardiac Infarction 

15 These examples were carried out on two 55-year-old male subjects with an 

apparatus and procedure as described in Example 1 above. Figure 6 illustrates quasi- 
stationary QT-RR interval hystereses for the male subjects. The curves fitted to the 
squares and empty circles relate to the first individual and illustrate a case of cardiac 
ischemia also detectable by the conventional ECG - ST segment depression 

20 technique. The curves fitted to the diamonds and triangles relate to the other subject, 
an individual who previously had experienced a myocardial infarction. These subjects 
exercised on a treadmill according to a quasi-stationary 20-minute protocol with a 
gradually increasing and gradually decreasing exercise load. Fitting curves are 
obtained using third-order polynomial functions. These cases demonstrate that the 

25 method of the present invention allows one to resolve and quantitatively characterize 

i 

the difference between (1) levels of ischemia that can be detected by the conventional 
ST depression method, and (2) low levels of ischemia (illustrated in Figure 5) that are 
subthreshold for the conventional method and therefore undetectable by it. The levels 
of exercise-induced ischemia reported in Figure 5 are significantly lower than those 
30 shown in Figure 6. This fact illustrates insufficient resolution, of a conventional ST 
depression method in comparison with the method of the present invention. 

/ 

% 
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EXAMPLE 6 
Hysteresis Curves in the Same Subject 
Before and After a Regular Exercise Regimen 
This example was carried out with an apparatus and procedure as described in 
5 Example 1 above. Figure 7 provides examples of quasi-stationary hystereses for a 55 
year-old male subject before and after he engaged in a practice of regular aerobic 

* 

exercise. Both experiments were performed according to the same quasi-stationary 
20-minute protocol with a gradually increasing and gradually decreasing exercise 
load. Fitting curves are obtained using third-order polynomial functions. The first test 
10 shows a pronounced exercise-induced cardiac ischemic event developed near the peak 
level of exercise load, detected by both the method of the present invention and a 

» 

conventional ECG-ST depression method. The maximum heart rate reached during 
the first test (before a regular exercise regimen was undertaken) was equal to 146. 
After a course of regular exercise the subject improved his cardiovascular health, 

15 which can be conventionally roughly, qualitatively, estimated by a comparison of 
peak heart rates. Indeed, the maximum heart rate at the peak exercise load from the 
first experiment to the second decreased by 16.4%, declining from 146 to 122. A 
conventional ST segment method also indicates the absence of ST depression, but did 
not provide any quantification of such an improvement since this ischemic range is 

20 sub-threshold for the method. Unlike such a conventional method, the method of the 
present invention did provide such quantification. Applying the current invention, the 
curves in Figure 7 developed from the second experiment show that the area of a 
quasi-stationary, QT-RR interval hysteresis decreased significantly from the first 
experiment, and such hysteresis loop indicated that some level of exercise-induced 

25 ischemia still remained. A change in the shape of the observed composite dispersion- 

restitution curves also indicates an improvement since it changed from a flatter curve, 

* 

similar to the flatter curves (with a lower excitability and a higher threshold, v r =0.3 to 
0.35) in Figure 2, to a healthier (less ischemic) more convex-shape curve, which is 
similar to the lower threshold curves (v r =0.2 to 0.25) in Figure 2. Thus, Figure 7 
30 demonstrates that, due to its high sensitivity and high resolution, the method of the 
present invention can be used in the assessment of delicate alterations in levels of 



■ 

» 
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cardiac ischemia, indicating changes of cardiovascular health when treated by a 
conventional cardiovascular intervention. 

a 

EXAMPLE 7 

r 

5 Calculation of a Quantitative 

Indicium of Cardiovascular Health 

This example was carried out with the data obtained in Examples 2-6 above. 

* 

Figure 8 illustrates a comparative cardiovascular health analysis based on ischemia 
assessment by the method of the present invention. In this example an indicium of 
10 cardiovascular health (here designated the cardiac ischemia index and abbreviated 

"CII") was designed, which was defined as a quasi-stationary QT-RR interval 

. ... 

hysteresis loop area, 5, normalized by dividing it by the product (rRR.max- 
^RR,min)(3QT,max-7QT,min)- F°r each particular subject this factor corrects the area for 
individual differences in the actual ranges of QT and RR intervals occurring during 

15 the tests under the quasi-stationary treadmill exercise protocol. We detennined 
minimum and maximum CII in a sample of fourteen exercise tests and derived a 
normalized index <CII> = (CII - CII ni i I1 )/(CII m ax - Cn,™) varying from 0 to 1. 
Alterations of <CH> in different subjects show that the method of the present 
invention allows one to resolve and quantitatively characterize different levels of 

20 cardiac and cardiovascular health in a region in which the conventional ST depression 
method is sub-threshold and is unable to detect any exercise-induced ischemia. Thus, 
unlike a rough conventional ST-segment depression ischemic evaluation, the method 
of the present invention offers much more accurate assessing and monitoring of small 
variations of cardiac ischemia and associated changes of cardiac or cardiovascular 

■ 

25 health. 

■ 

* m 

i 

EXAMPLE 8 

■ 

Illustration of Rapid Sympatho-Adrenal Transients 
Figure 9 illustrates a typical rapid sympathetic/parasympathetic nervous and 
30 hormonal adjustment of the QT (panels A, C) and RR (panels B,D) intervals to an 
abrupt stop after 10 minutes of exercise with increasing exercise load. All panels, 
depict temporal variations of QT/RR intervals obtained from the right precordial lead 
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■ V3 of the 12-lead multi-lead electrocardiogram. A sampling rate with which a 

> 

waveform analyzer determined QT and RR intervals was equal to 15 samples per 
minute. A human subject (a 47 years-old male) was at rest the first 10 minutes and 

♦ 

then began to exercise with gradually (during 10 minutes) increasing exercise load 
5 (Panels A, B - to the left from the RR, QT minima). Then at the peak of the exercise 
load (heart rate about 120 beat/min) the subject stepped off the treadmill in order to 
initialize the fastest RR and QT interval's adaptation to a complete abrupt stop of the 
exercise load. He rested long enough (13 minutes) in order to insure that QT and RR 
intervals reached post-exercise average stationary values. Panels C and D demonstrate 
1 0 that the fastest rate of change of QT and RR intervals occurred immediately after the 
abrupt stop of the exercise load. These rates are about 0.015 s/min for QT intervals 

* 

while they vary from 0.28s to 0.295s and about 0.15s/min for RR intervals while they 

i 

grow from 0.45s to 0.6s. Based on the above-described experiment, a definition for 

■ 

"rapid sympathoadrenal and hormonal transients" or "rapid autonomic nervous 
1 5 system and hormonal transients" may be given. 

• • * 

Rapid transients due to autonomic nervous system and hormonal control refer 
to the transients with the rate of 0.1 5 s/min for RR intervals, which corresponds to the 
heart rate's rate of change of about 25 beat/min, and 0.02s/min for QT intervals or 
faster rates of change in RR/QT intervals in response to .a significant abrupt change 
20 (stop or increase) in exercise load (or other cardiac stimulus). The. significant abrupt 
changes in exercise load are defined here as the load variations which cause rapid 
variations in RR/QT intervals, comparable in size with the entire range from the 
exercise peak to the stationary average rest values. 

25 EXAMPLE 9 

Illustration of a Quasi-Stationary Exercise Protocol 
Figure 10 illustrates a typical slow (quasi-stationary) QT (panel A) and RR 
(panel B) interval adjustment measured during gradually increasing and gradually 
decreasing exercise load in a right pre-cordial V3 lead of the 12 lead 
30 electrocardiogram recording. The sampling was 15 QT and RR intervals per minute. 
A male subject exercised during two consecutive 10 minute long stages of gradually 
increasing and gradually decreasing exercise load. Both QT and RR intervals 
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* 

gradually approached the minimal values at about a peak exercise load (peak heart 
rate -120 beat/min) and then gradually returned to levels that were slightly lower than 
their initial pre-exercise rest values. The evolution of QT and RR intervals was well 
approximated by exponential fitting curves shown in gray in panels A and B. The 
5 ranges for the QT-RR interval, there- and-back, time variations were. 0.34s - 0.27s - 

* * 

0.33s (an average rate of change ~0.005s/min) and 0.79s - 0.47s - 0.67s (an average 
rate of change ~6.032s/min or ~6 beat/min) for QT and RR intervals, respectively. 

* 

The standard root-mean-square deviation, cr, of the observed QT and RR intervals, 
shown by black dots in both panels, from their exponential fits were on an order of 
10 magnitude smaller than the average difference between the corresponding peak and 

t 

m 

rest values during the entire test These deviations were cr~0.003s for QT and <r~0.03s 
for RR intervals, respectively. According to Figure 9 (panels C, D) such small 
perturbations, when associated with abrupt heart rate changes due to physiological 
fluctuations or due to discontinuity in an exercise load, may develop and decay faster 
15 than in 10s, the time that is 60 times shorter than the duration of one gradual 

■ 

(ascending or descending) stage of the exercise protocol. Such a significant difference 
between the amplitudes and time constants of the QT/RR interval gradual changes and 
. abrupt heart rate fluctuations allows one to average these fluctuations oyer time and 
fit the QT/RR protocol duration dynamics by an appropriate smooth exponential-like 
20 function with a high order of accuracy. A simultaneous fitting procedure (panels A, 

■ 

B) determines an algorithm of a parametrical time dependence elimination from both 
. measured QT/RR data sets and allows one to consider QT interval for each exercise 
stage as a monotonic function. 

Based on the above-described experiment a definition for a gradual, or "quasi- 
25 stationary" exercise (or stimulation) protocol, can be quantitatively specified: A 
quasi-stationary exercise (or stimulation) protocol refers to two contiguous stages 
(each stage 3, 5, 8 or 10 minutes or longer in duration) of gradually increasing and 

gradually decreasing exercise loads or stimulation, such as: 

■ 

1. Each stage's duration is approximately an order of magnitude (e.g., at least 
30 about ten times) longer than the average duration (~ 1 minute) of a heart rate 

adjustment during an abrupt stop of the exercise between average peak load 



*■ • 



» 
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rate (~ 120 - 150 beat/min) and average rest (~ 50 - 70 beat/rain) heart rate 
values. 

< t m 

2. The standard root-mean-square deviations of the original QT/RR interval 

• a 

data set from their smooth and monotonic (for each stage) fits are of an order 
5 of magnitude (e.g., at least about ten times) smaller than the average 

differences between peak and rest QT/RR interval values measured during the 
entire exercise under the quasi-stationary protocol. 
As shown above (Figure 10) a gradual quasi-stationary protocol itself allows one to 

* 

substantially eliminate abrupt time dependent fluctuations from measured QT/ RR 
10 interval data sets because these fluctuations have short durations and* small 

■ 

amplitudes. Their effect can be even further reduced by fitting each RR/QT interval 
data set corresponding to each stage with a monotonic function of time. As a result the 
fitted QT interval values during each exercise stage can be presented as a substantially 
monotonic and smooth function. of the quasi-stationary varying RR interval value. 

15 Presented on the (RR-interval, QT-interval)-plane this function gives rise to a loop, 
whose shape, area and other measures depend only weakly on the details of the quasi- 
stationary protocol, and is quite similar to the hysteresis loop presented in Figure 2. 
Similar to a generic hysteresis loop, this loop can be considered as primarily 
representing electrical conduction properties of cardiac muscle. 

20 It is well known that exercise-induced ischemia alters conditions for cardiac 

electrical conduction. If a particular individual has an exercise-induced ischemic 
event, then one can expect that the two experimental composite dispersion-restitution 
curves corresponding to the ascending and descending stages of the quasi-stationary 
protocol will be different and will form a specific quasi -stationary hysteresis loop. 

w * 

* < 

25 Since according to a quasi-stationary protocol the evolution of the average values of 
QT and RR intervals occurs quite slowly as compared with the rate of the transients 
due to sympathetic/parasympathetic and hormonal control, the hysteresis loop 
practically does not depend on the peculiarities of the transients. In that case such a 
hysteresis can provide an excellent measure of gradual ischemic exercise dependent 

* 

30 changes in cardiac electrical conduction and can reflect cardiac health itself and 
cardiovascular system health in general. 
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It should be particularly emphasized that neither J. Sarma et al, supra, nor A. 
Krahn's et al. supra, report work based on collecting quasi-stationary dependences, 
which are similar to the QT- interval -RR interval dependence, since in fact both 
studies were designed for different purposes. To the contrary, they were intentionally 

* ■ * * 

5 based on non - quasi-stationary exercise protocols that contained an abrupt exercise 
stop at or near the peak of the exercise load. These protocols generated a different 
type of QT/RR interval hysteresis loop with a substantial presence of non-stationary 
sympatho-adrenal transients (see Figure 9 above). Thus these prior art examples did 
not and could not include data which would characterize gradual changes in the 
10 dispersion and restitution properties of cardiac electrical conduction, and therefore 
included no substantial indication that could be attributed to exercise-induced 
ischemia. 

■ * 

Examples 10-12 

15 , . Data Processing at Steps 44-45a,b (Figures 4Aand4B) 

The following Examples describe various specific embodiments for carrying out the 
processing, shown in Figure 4B, where the software implemented at steps 44b and 45b 
performs the following major steps: 

■ 

(i) Generates sufficiently smooth time dependent QT and RR data sets by 
20 averaging/filtering and fitting these averaged data by exponential or any other 

* 

• * 

suitable functions; 

(ii) Combines into pairs the points of RR and QT data sets that correspond to the 
same time instants, thereby generating smooth QT/RR curves for the ascending and 
descending heart rate stages (branches) on the (rRK,7QT)-plane or a similar plane or 

* 

25 its image in computer memory; 

(Hi) Closes the ends of the QT/RR curves to transform them into a closed loop 
and determines an area, S, or a similar measure of the domain confined by the said 
loop and computes an index to provide a quantitative characteristic of cardiovascular 

■ 

health in a human subject based on a measure of this domain. 
30 . Each of these three major steps may consist of several sub steps, as illustrated in each of 
the methods shown in Figure 11, Figure 13 and Figure 15, Figure 16 and as discussed in 
greater detail below. 
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Example 10 

A Method of Optimized Consolidation of a Moving Average, 

Exponential and Polynomial Fitting 
5 1. Raw data averaging/filtering (box 1 in Figure 11). The raw data set 

consists of two subsets {t^T^} i=l,2,.. JVrr-1, and {t l QT ,TQ T } 9 /=1,2,.. JVqt-1, 

where t and' T x denote the i-th sampling time instant and the respective, RR or QT, 

interval duration (subsript x stands for RR or QT). In order to simplify the notation we 
shall omit the subscripts RR and QT when it is applicable to both sub-samples. It is 

10 convenient to represent each data point as a two-component vector u\=(t\T). The 
filtering procedure in this example comprises a moving averaging of neighboring data 
points. We shall denote a moving overage over a set of adjacent points by angular 
brackets with a subscript indicating the number of points included in the averaging 
operator. Thus, the preliminary data filtering at this sub step is described by the 

15 equation 

■ 

Wa-^-y- 1 ^." 001) 

for each /=1,2,.. JV-1, where N is a number of data points in the corresponding (RR or 
QT) raw data set. Thus, the RR and QT interval durations and the corresponding 

■ 

20 sampling time points are identically averaged in order to preserve the one-to-one 

correspondence between them. This procedure removes the high frequency noise 

* 

present in the raw data. This preliminary smoothing can also be described in the 
frequency domain and alternatively achieved via appropriate low-pass filtering. All the 
following processing pertinent to this example will be done on the averaged data points - 
25 <Wi>2 and the angular brackets will be omitted to simplify the notation. 

2. Preliminary estimation of the t-coordinates of the minima (box 2 in Figure 

11), The sub step consists of preliminary estimation of the time instants, t*£ R and fjjj.of 

the minima of the initially averaged RR-interval and QT-interval-data sets, 
respectively. The algorithm does a sequential sorting of the data sets choosing M data 
30 points corresponding to M least values of the RR- or QT-intervals and then averages the 
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result. (In our examples, M=10.) First, the algorithm finds the u-vector w. =(/ / *,r /| ) 

corresponding to the shortest interval T h . Next, this data point is removed from the. 
data set {u\} and the algorithm sorts the remaining data set and again finds a u-vector 

u i2 = (t' 2 ,T' 2 ) corresponding to the shortest interval T 2 . The minimum point is again 

* * 

5 removed from the data sets and sorting is repeated over and over again, until the A/* u- 

• * - 

vector u i = (t lm ,T' a ) corresponding to the shortest interval T' m is found. Next, the 

algorithm averages all M pairs and calculates the average minimum time coordinate t 
defined as 

* * + 

M 

io . . (10-2) 

> 

Finally, the algorithm determines the sampling time instant t im , which is nearest to t. 
We thus arrive to t 1 ^ and t J £ T for RR and QT data sets, respectively. 

* 

3. The first correction of the coordinates of the preliminary minima (box 3 in 
Figure 11). At this sub step the first correction to the minimum coordinates and 

1 5 t ] g r is found. This part of the algorithm is based on the iterative exponential fitting of 

■ 

the 7-component of the initially filtered data set {ui]=f{t \T*} by functions of the form 

< 

T(t) = Aexj>US\t-t m \], (10.3) 

t 

20 where / m is the time instant of the corresponding minimum determined at the previous 
sub step {i.e., t m for RR-intervals and t m =t J 0 m T for QT-intervals). The fitting is 

m 

i 

done separately for the descending (t<f) and ascending (t>t m ) branches of u(t) with the 
same value of constants A and t m and different values of P for both branches. The initial 

value of A is taken from the preliminary estimation at sub step 2: A=T in and 

25 t m = & {Le., A^T^,t m = fj- for RR-intervals and A = T£ ,t m = t J g T for QT- 

intervals). The initial value of constant p. for each branch and each data sub set is 

■ 

obtained from the requirement that the initial mean root square deviation cr=a 0 of the 

* 

entire corresponding branch from its fit by equation (10.3) is minimum. In each 
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subsequent iteration cycle new values of constants A and P are determined for each 
branch. At the beginning of each iteration cycle a constant A is taken from the previous 
step and the value of p is numerically adjusted to minimize the deviation value, o. Such 
iterations are repeated while the mean square deviation o is becoming smaller and are 
5 stopped when a reaches a minimum value, 0=0^. At the end of this sub step the 
algorithm outputs the corrected values A and p and the respective values of cr££ and 
■ a QT 

u min • 

4. The preliminary smooth curves (box 4 in Figure 11). The sub step consists 
of calculating a series of moving averages over p consecutive points for each data 
10 subset {w} as follows 



P M 



(10.4) 



We shall refer to the quantity p as the width of the averaging window. The calculation 
15 is performed for different values of p, and then the optimum value of the width of the 
averaging window p=m is detennined by minimizing the mean square deviation of the 

■ 

set {U?} = {t^T*} from its fit by Eq. (10.3) with the values of parameters A, F, and fi 
determined at sub step 3. Having performed this procedure for each component of the 

■ 

data set (RR and QT) we arrive to the preliminary smoothed data set {uf} = {*",2 , / "} , 
20 where 



i+m-l 



(10.5) 



The number of points in such a data set is N m =N-m+l, where N is the number of data 

■ 

25 points after the initial filtering at sub step 1 . 

5. Correction of the preliminary smooth curves and the second correction of 
the QT and RR minimal values and its coordinates (box 5 in Figure 11). At this sub 
step we redefine the moving averages (using a smaller averaging window) in the 
vicinity of the minimum of quantity T (RR or QT interval). This is useful to avoid 
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distortions of the sought (fitting) curve T=T(t) near its minimum. The algorithm first 

. i 

■ • m 

specifies all data points {t,T 1 } such that T x lie between minir 1 } and minfr'j-f-c^ . 

These data points have the superscript values f from the following set 

/ = {/:7U[mm{r.} J min{r}+a min ]}. (10.6) 

5 Let us denote by i 0 e / and i e / the subscript values corresponding to the earliest and 

* » 

i 

latest time instants among {t\} (iel). We thus set i<f=mm{I} and z q =max{7}. Denoting 
the number of points in such a set by #rr for the RR data and by #qt, for the QT data 
we determine the width, q, of the averaging window as the minimum of the two, 
g=min{gQT,<?RR} . Now we write the moving average for the data points in the vicinity 
10 of the minimum as 

■ 

* 

q h 

* * 

Such modified averaging is applied to all consecutive data points u x with the subscript i 
15 ranging from z 0 to i g . The final (smoothed) data set consists of the first z"o-l data points 

defined by Eq.(10.4) with i = l,2,...i 0 -1, q points defined by Eq.(10.7) with 

* = Wo + ™d N-i q -m points defined again by Eq.(10.4) with 

i-i g +l,i q + 2,...,iV-7H + l. This the final set {u. } can be presented as 

20 {u.} = {uf :i6{l,2,...,/ 0 -l}U{i, +l,../-)«+l}}U{« ( ' :i'e0 o .'o + L -,',}> (10-8) . 

■ 

At the end of sub step 1 6.5 the algorithm detennines the final minimum values of QT 
and RR interval and the corresponding time instants. The algorithm sorts the smoothed 

■ * 

data sets (10.8) and determines u ) corresponding to the minimum value 

25 of Tt . This can be written as 

4 

* 

«— s {0:,Ti.)-TL = nWf *}} = *<,. (10.9) 



V 
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* 

6. Final fitting and the final smooth QT and RR curves (box 6 in figure 11). 

At this sub step the parameters of functions representing final smooth QT and RR 

curves are found. First, each data set {«/} is split into two subsets {u* } and {u\} 

5 corresponding to the descending (/,■ < ti m ) and ascending (ti>t fm ) branch, respectively. 

We also shift the time origin to the minimum point (i^ ) and change the axis direction 
of the descending branch. These redefined variables are thus given by : 



10 



{u'} = {(t i ,f) S (t, m -t i ,T i ):iZi m }}, 



{ul} = {(tlj;) = (t i -t la ,T l ):i*i m }. 



(10.10a) 
(10.10b) 



15 



« 

Notice that the minimum point is included in both branches. Next we perform a linear 

i 

regression by fitting each branch with a 4-th order polynomial function in the following 
form 

T ± (t ± ) = a ± t 4 ± +bji +c ± yl +d (10.11) 



The linear term is not included into this expression since this function must have a 
minimum at 4=0..The value of is defined by 



20 



d = T 



nun 



. (10.12) 



25 



For computational convenience the variables w± are transformed into z± — T±-d and 
then the coefficients a±, &±, and c± are determined from the condition that expression for 



the error 



£ =I>i -a±('i) 4 -cMff ' 



(10.13) 



where summation is performed over all points of the corresponding branch, is 
minimized. We thus obtain four similar sets of linear algebraic equations - two for 

* 

each of the u ± branches of the two RR and QT data sets. These equations have the 

♦ 

form 
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lei) 8 


lei) 7 


lei) 6 






Eei) 4 *i 




lei) 6 


lei) 5 








lei) 6 


lei) 5 


leir 






Iei) 2 *i 



(10.14) 



10 



15 



20 



25 



All summations are performed over all points of the branch. Now the QT and RR 
interval time dependences curves can be presented as the follows 

» 

and 



min 



(10.15) 



QT 

r + (t-t QT ) t>t Q1 



QT 
min 



(10.16) 



where 7^(?) and Tg T {t) are the fourth order polynomials given by Eq.(lO.ll) and 



and fj? are the time coordinates of the corresponding minima defined by Eq. (10.9). 

Thus, formulas (10.10)-(10.16) determine the final smooth QT and RR curves with 
minima defined by formula (10.9). 

7. Final smooth hysteresis loop (box 7 in Figure 11). At this sub-step the 

* 

software first generates a dense (N+l)-point time-grid x k = t siui + k(t m & - ts^/N, where 
A=0,l,2,. ..^V (AT=1000 in this example) and t stzn and / e nd are the actual values of time at 
the beginning and end of the measurements, respectively. Then it computes the values 

of the four functions 7]£ (l - f ) and 7^ r (x - f^T ) on the grid, which parametrically 



represents the final smooth curves ( (x - ) , Tq T (x — ^ ) ) in computer memory 

for the ascending ( + ) and descending ( - ) branches, respectively. This procedure is a 
computational equivalent of the analytical elimination of time. Next, the software also - 

* » 

plots these smooth curves on the (7RR,7QT)-plane or on another, similar plane, such as 
(/rr*'7qt), where /=1/7rr is the instantaneous heart rate. Finally, the algorithm adds to 
the curves a set of points representing a closing line, which connects the end point of 
the lower (descending) branch with the initial point on the of the upper (ascending) 
branch and, thus, generates a closed QT/RR hysteresis loop. 
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8. A measure of the domain bounded by the QT/RR hysteresis loop (box 8 in 
Figure 11). At this sub step a measure of the "domain inside the QT/RR hysteresis loop 
is computed by numerically evaluating the following integral (see definition above): 

5 • S^jlpiT^T^dT^dTgj. (10.17) 

where Q. is the domain on the (rRR,rQi)-plane with the boundary formed by the closed 
hysteresis loop, and p(7rr,7qt) is a nonnegative (weight) function. In this example we 
can take p(7rr,7qt) =1 so that S coincides with the area of domain D, ,or we set p = 
1/(7rr) so S coincides with the dimensionless area of the domain inside the hysteresis 
1 0 loop on the (/]TQT)-plane, where /=1/Trr is the heart rate. 

' 9 

9 ■ 

I 

■ 

EXAMPLE 11 
A method of a sequential moving average 
15 1. Raw data averaging/filtering* The sub step (similar to 10.1) consists of raw 

■ 

data averaging (filtering) according to formula (10.4). This is performed for a 
preliminary set of values of the averaging window width, p. 

2. Final moving average smoothing^ The sub step includes subsequent final 

smoothing of the preliminary smoothed data represented by the set {uf } found at the 
20 previous step: as follows. 

U'mluf) =— (11.1) . 

* m f j=i 

The window width mf is varied numerically to achieve optimum smoothness and 
accuracy of the fit. The optimally averaged data points form smooth curves on the 

« 

corresponding planes (Figure 12,14). 
25 3. A measure of the domain inside the QT/RR hysteresis loop. The procedure 

on this sub step is as defined in Example 1 0, step 8, above. 

■ 

* * 

* 

* 

EXAMPLE 12 

30 A Method of Optimized Nonlinear Transformations 
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* 

Figure 15 illustrates major steps in the data processing procedure involving 
our nonlinear transformations method. The first three stages for the RR and QT data 
sets are quite similar and each results in the computation of the fitted bird-like curves, 
7 rr = r RR^ md r QT = jQT(^ m a dense time-grid. Having computed both bird-like 

5 curves one actually completes the data processing procedure, because after 
appropriate synchronization, these two dependences parametrically along with a 

closing line represent the sought hysteresis loop on the (t** , T QT ), or similar, plane. 

0 

We will describe our method in general terms equally applicable for the QT and RR 
interval data sets, and indicate the specific instants where there is a difference in the 
10 algorithm. 

1. Preliminary data processing stage. Let {7k}, £=1,2, "Jf, be a set of the 
measured RR or QT interval durations, and be a set of the corresponding time 

instants, so that t\ and fa are the starting and ending time instants (segments) of the 

■ 

entire record. The first three similar stages (stages 1 through 3 in Figure 15) are the 
15 preliminary stage, the secondary, nonlinear-transformation stage, and the 
computational stage. A more detailed data flow chart for these stages is shown in 
Figure 16. The preliminary stage, indicated by boxes 1 through 7 in Figure 16, is a 
combination of traditional data processing methods and includes: smoothing 

* * 

(averaging) the data sets (1), determination of a region near the minimum (2), and 
20 fitting a quadratic parabola to the data in this region (3), checking consistency of the 
result (4), renormalizing and centering the data at the minimum (5), cutting off the 
data segments outside the exercise region' and separating the ascending and 
descending branches (6), and finally filtering off a data segment in near the minimum 
(7). Let us discuss these steps in more detail. 
25 Box 1 in Figure 16, indicates the moving averaging and itself consists of two 

steps. First, we specify the initial value of the moving average procedure parameter m, 

« 

* - m 

as - 



m = max< r 



[ 



nun 



(12.1) 



« 



WO 02/00113 PCTAJS01/20306 



25 



-45- 

where r{x) is the integer closest to x, and the values of N c and m m \ n are chosen 
depending on the number N of data points and the amount of random fluctuations in 
the data (the size of the random component in the data). In our examples we have 

chosen iV c =100 and m m i n =3. The value of m can be redefined itefatively at the later 

, * » 

stages and its choice will be discussed therein. Next we compute the moving average 
for ft} and {T{} data sets with the given averaging parameter m as follows: 

m w m , =l 



The subscript m will be omitted below if m is fixed and no ambiguity can arise. The 
1 0 next step in our algorithm is the initial determination of a time interval on which the 
parabolic fit will be performed, (the parabolic fit interval). This step is represented in 
Figure 16 by Box 2. We note that this data subset can be redefined at a later stage if a 
certain condition is not satisfied. In the current realization of our algorithm this 

region defined differently for faTf*} and {t i9 T? T } data sets. For the data set 

* 

15 {tuTf* } the initial parabolic fit segment is defined as'the data segment {fcrf* } with 

■ 

t 

all sequential values of i between i\ and h where U=n m \ n -Am and / 2 = "min + Am, and the 
integer parameters n m { n and Am are defined as follows. The number rt m in determines 

• » 

. the time instant t n among the original set of time instants which is the closest to the 
minimum on the averaged data set {<iY>,<7j>} . In other words, t is the nearest to 

■ 

20. the average time instant <^m> which corresponds to the minimum value of the average 
RR-intefval <Ti RR >, that is, the time instant with the subscript value M defined by the 

* 

condition 

* 

4 

<T£ R >= m m{<T ! «* >} - (12.3) 



The value of Am is linked to the value of m by the condition A/n=2m. The link 
between Art and m arises from the requirement that the algorithm is stable and 
consistent. The consistency is the requirement that the positions of the minimum of 
the curve obtained by moving averaging and by quadratic fitting were approximately 
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9 

the same. On the other hand, it is also important that the value of m is sufficiently 

• • • 

small because quadratic fitting is done with the purpose to describe the. data only 
locally in the immediate vicinity of the minimum. 

For the data set {t\,Tx} the initial parabolic fit segment is defined as the data 
subset that consists of all consecutive points belonging to the lower portion of the 
non-averaged {T, QT } data set. We thus define i\ as the first (minimum) number i such 

■ 

♦ 

that 



T, ST < min (rf )+ R 



max (rf )- min (rf ) 



(12.4) 



where R is a parameter, 0<i?<l, that determines the portion of the data to be fitted 
with the quadratic parabola. In our calculations we set £=1/8=0.125, so we fit with the 
parabola the data points that correspond to the bottom 12.5% of the interval durations. 
Thus, the subscript z'i is the first (smallest) value of i such that condition (12.4) is 
15 . satisfied. Similarly, i 2 is the last (greatest) value of i such that condition (12.4) is 
satisfied. The initial fit region is then defined as the following non-averaged data 

subset: { (/. , Tf r ):/ = /,, i x + 1, /, + 2, . . . , i 2 } . This method can also be used for 

determining an initial data segment for the quadratic polynomial fitting of the RR-data 
set. . • 

20 At the next step, shown by box 3 in Figure 12.2, we fit the data in this region 

(for each data set) by a parabola so that the data are approximately represented by the 
equation 



T k *PA+P 2 t k +Pi, k = J\,j\+l,j\+2,~'>U (12.5) 



This is done using usual linear regression on this data subset. At the next step (Box 4 
in Figure 16) we check if the parabola has a minimum (i.e., if P\>0). If this condition 
is satisfied, the determination of the parabolic fit interval and the fitting parabola itself 
is completed. Otherwise, we enter the loop indicated by boxes la through 3a in Figure 
30 16. The first step there is similar to the above second step utilized for the RR interval 
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data set and defines an extended segment for the parabolic fitting. This is done by 
replacing m with m+\ and using this new value of m to redefine Am -2m and also 
calculate new averaged data sets in accordance with Eq.(12.2) as indicated by Box la 
in Figure 16. This allows us to evaluate a new value of n m \ n and new values of /i=nmin- 
5 Am and / 2 = "min+Am with new value of m (Box 2a). then the parabolic fit is done 
again (Box 3a) and the condition P\>0, that the parabola has a minimum, is checked 
again. If it is satisfied, the determination of the fit interval and the quadratic fit 
coefficients procedure is completed. Otherwise we replace m with m+l and repeat the 
process again and again until the condition Pi>0 is satisfied. This condition ensures 

10 that the extremum of the quadratic parabola is a minimum indeed. By the very design 
of the exercise protocol the average heart rate reaches a certain maximum somewhere 
inside the load stage . and therefore corresponding average RR-interval reaches a 
minimum. This ensures that a so-defined data segment exists and is unique, and 
therefore the coefficients of the quadratic parabola are well defined. In short, we use 

15 the shortest data segment that is centered around the minimum of the averaged 
data set and that generates a quadratic parabola with a minimum (Pi>0). . 

The parabolic fit defines two important parameters of the data processing 

> • * * 

procedure, the position (AnimTmin) °f th e minimum on the (/,J)-plane as. follows 



20 ^ = "^' ^^"Ir (12 ' 6) 



These parameters are final in the sense that our final fit curve will always pass 
through and have a minimum at the point (fmim^min)- Th e coordinates (fmin^min) thus 
constitute parameters of our final fit bird-like curve. Having found the ordinate T m i n of 
25 the minimum we can renormalization of the T-data set as follows: 



02.7) 



min 



We also take the abscissa of the parabola's minimum t m m as the time origin and define 
30 the time components of the data points as follows 
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h = - Ctm > *i 



These two data transformations are indicated by Box 5 in Figure 16. We also restrict 
5 the data set only to the exercise period, plus some short preceding and following 
intervals (Box 5). The conditions f7 ^ 0, andt*j > o define the descending, {tJ,T\ }, 

and ascending, {t*j, 7]}, branches, respectively, so the corresponding data points can 

be readily identified and separated (Box 6). Given the durations and f a of the 
descending and ascending load stage, respectively, we can reduce the original set by 

10 cutting off the points on the descending branch with, tT < - r d and the points on the 

* 

ascending branch with tt > f a (Box 6 in Figure 16). This determines the minimum 
value io of subscript i for the descending branch, and the maximum value y max of 
subscript j on the ascending branch. 

The final step of the preliminary data processing is the conditional sorting 

* 

15 (Box 7). The conditional sorting removes all consecutive points such that at least one 
of them falls below the minimum of the parabola. The removal of the points below the 
minimum is necessary because the nonlinear transformation is possible only when 

* * * - 

yi>l . Eliminating only separate points below the minimum y-1 would create a bias in 
the data.. Therefore, we have to eliminate an entire segment of the data in the vicinity 
20 of the minimum. It should be remembered though, that these points have already been 
taken into account in the above quadratic filtering procedure that determines (ton Jmin) 
via Eq, (12.6). " 

Thus, the preliminary data processing results in two data sets corresponding to 
the descending and ascending load stages. The descending data set is defined as a set 

25 of consecutive pairs (tj , y$=( tj,yj) with i = i 0 , fo+1, *o+2, z max , where is 

determined by the conditional sorting as the largest subscript i value still satisfying 
the conditional. Similarly, the ascending data set is defined as a set of consecutive 
pairs (fiVi) = faVi*) withy = y m in, 7o + l , jo+2, "Wmax, where y min is determined by the 
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10 



conditional porting as the first subscript j value on the ascending branch starting- from 
which the condition^ >1 is satisfied for all subsequent jr. 

2. Secondary data processing. At the second stage we introduce a 
fundamentally new method of nonlinear regression by means of two consecutive 
optimal nonlinear transformations. The idea of the method is to introduce for each 
branch two appropriate nonlinear transformations of the dependent and independent 
variables, y=flu) and w=Hp(0, both transformations depending on some parameters, 
and choose the values of the parameters in such a way that the composition of these 
transformations J[(p(t\)) would provide an approximation for y\. Let (t\ ' , yf) and 
yt) be the cut, conditionally sorted and normalized data set corresponding to the 
descending and ascending branches representing the dynamics of the RR- or QT- 
intervals during exercise. The nonlinear transformation u is defined via a smooth 
function 



15 



y=fy (")> 



(12.9) 



that has a unit minimum at u = 1,'so that f y (l) = 1, //(1) = 0, and f y (u) grows 

monotonously when u > 1 and decreases monotonously when u < 0. The subscript y 
represents a set of discrete or continuous parameters and indicates a particular choice 

* • * 

■ • 

20 of such a function. Let us denote by f~ (w) the monotonously decreasing branch of f r 
(u) and its monotonously increasing branch, by / Y + (u) . We can thus write 



/,(«)= 



f~(u),when u<l, 
f* \u\when u > 1. 



(12.10) 



» 

25 Let u = g* (y) and u = g y " (y) be the inverse functions for the respective branches of 
f y (u) . The functions g y + (y) and f* (u) are monotonously increasing, while the 

■ 

* < ■ 4 

■ functions g~(y) and f~(u) are monotonously decreasing. Let {t., yj} represent the 
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> 

data for the descending segment of the data set, i.e. t. < , and {t., yj t } represent 
the data for the ascending segment of the data set i.e. for t. > t. . The transformation 

w <*j i nun 

* 

?;,=g;(y,-) . (i2.il) 

5 

* maps the monotonously decreasing (on the average) data set {t n y7} into a 
monotonously increasing one: 

« 

^y;}^^g;{y?h^;A ■ - 02.12) . 

■ B * 

1 

10 

■ 

Moreover, the average slope of the original data set is decreasing as t\ approaches / m m 

■ 

and eventually vanishing at the minimum. In contrast, the average slope of the 

- 

transformed data set is always nonzero at /=/ m i n . Similarly, the transformation 
15 <j=gtW (12.13) 

* 

1 

maps the monotonously increasing (on the average) data set {t., yj) into a 
monotonously increasing one:- 

* 

In our examples we used a discrete parameter y taking two values y -1 and 7 =2 that 
correspond to two particular choices of the nonlinear y-transformation. The first case 
(Y=l) is described by the equations 

25 

y=7iW-l + («-iy: ^=g;(y)^y-^l 9 u\ = gy{y) = y + ft^- 
(12.15) 

The second case, y=2, is described- by the equations 
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y /,(«)■«+-; u=gl(y)^\-^A, u*=gl{y) = \ + &T. (12.16) 



• » 



Both functions reach a minimum y=l at w=l . An example of such a function 
5 . depending on a continuous parameter y >0 is given by 



/rM a ^TT" + 57~7 K + r + (12.17) 



The parameter 6 Y has been determined from the condition that/y(w) has a minimum at 

■ 

1 0 w=l and the coefficient A Y is determined by the condition ^(1)=1, which yields 



7 ~~ ..3 , „2 ~~ 



+ y +2y 



(12.18) 



In our numerical examples below we utilize the discrete parameter case with y taking 
15 two values, 1 and 2. The original and transformed sets are shown in Figures 17 (RR 
intervals) and 18 (QT-intervals). Panels A in both figures show the original data sets 
on the (r,7)-plane and panels B and C show the transformed sets on the (r,w)-plane, for 
y - 1 and for y = 2, respectively. The parabola minima on panels A and C are marked 
with a, circled asterisk. The original data points concentrate near a non-monotonous 
20 curve (the average curve). The data points on the transformed plane concentrate 
around a monotonously growing (average) curve. Moreover, the figure illustrates that 
the transformation changes the slope of the average curve in the vicinity of t = f min 
from zero to a finite, nonzero value. On the (M m i n ,w)-plane the point corresponding to 
the parabola's minimum is (0,1). This* point is also marked by a circled asterisk on 
25 panels' B and.-C of Figures 17 and 18. 

+ 

Let us introduce a pair of new time variables x" and t for the descending and 
ascending branches, respectively. We shall count them off from the abscissa of the 

p ■ 

minimum of the fitted parabola and set 
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T J =t j ~ 'nun ' 0 -'min»7 =1 ' 2 >"-.-/ + 



where /" and J + are the number of data points on the descending and ascending 
branches, respectively. The (formerly) descending branch can be treated in exactly the 
5. same way as the ascending one if simultaneously with the time inversion given by the 
first line in Eq. (12.19) we perform an additional transformation of the descending 
branch ordinates as follows 



v k =!-«;. (12.20) 



In these variables both branches (x[*,v,) and (x/^Wj*) are monotonously growing and 

start from the same point (x=0,w=l) at which both have nonzero slope and possess 
similar behavior (convexity). 

Since both branches are treated in exactly the same way, we shall simplify the 
15 notation and temporarily omit the superscripts ± and write (xk,w k ) for any of the pairs. 

(xr,Vi) or (xj*,^/*). We shall fit the data set {u k } with {<p(a , P , x k )}, that is represent 

{w k } as 



u k *<p{a,/3,K,x k ) (12.21) 



where the function <p depends linearly on K and is defined as 

<p{a,p,K,T) = \+KZ(a,$,t) (12.22) 

* ■ 

where 
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ai tY 
— rl+— 



s s a In 



a 



r 



1 + — 



1+ 



a 



-1 



when /?>0 
when ft = 0 



-l<p<0 



A family of functions ^ (oc,P,t) is shown in Figure 19 for fifteen values of parameter 

p. The parameter a is completely scaled out by plotting the function on the plane 
5 (x/a 3 £/a). The function £ (a,P,x) is continuous in all three variables and as a function 

of x .at fixed a and p has a unit slope at t=0, £ '(a,P,0)=l. Therefore, at t-»0 the 

function £ possesses a very simple behavior, £ ~r, which is independent of a and p 

(the size of the region of such behavior of course depends on a and P). The function 
i; possesses the following important feature: when P passes through the point P=0, its 

» ■ 

10 asymptotic behavior at t-*x> continuously changes from the power function, £ ~r p at 
P>0, through 4 ~ln(t) at p =0, to a power of the logarithm, % ~ln 1+p (x) when -1<P<0. 



15 



When p^-1, the behavior of § changes once more, and becomes % ~ln(ln(t)). The 

convexity of any function of the family is the same when P<1 . 

We shall explicitly express parameter K via a and P using the requirement for 
a given pair of a and P Eq-s. (12.22) and (12.23) ensure the best data fit in the u-space 
in a certain vicinity of the point (t=0,u=1). Let K\ be the number of points where the 
fit is required. The corresponding quadratic error is then a function of K, a and p that 
can be written as 



20 



* 

z< u \K,a,fl).= 'Z[u l -K<p{a,P,T k )] 

kSJC, 



(12.24) 



and the requirement of its minimum immediately yields the following expression for 
K 
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K-K{q,fi)m 



kiK, 



kZKt 



(12.25) 



. In our calculations we used such K\ value that would include all adjacent points with 
the values of u between u=l and u = 1 + 0.4 (u^ - 1). We have thus reduced the 
5 number of fitting parameters in our fitting procedure to two continuous parameters, a 

■ ► 

and P, and one discrete parameter, y. The fitting function thus becomes 



(12.26) 



10 The values of parameters a, P and K (and. y) are now directly determined by the 
condition that the fit error in the y-space . 



*J"W)-Zk -ft & +*MH*A* Jf 



(12:27) 



15 is niinimum. The sum in Eq. (12.27) is evaluated numerically on a grid (a, P ) values 

for y equal 1 and 2. Then the values of a, [5 and y that deliver a minimum to ej y) are 

found via numerical trials. Calculations can then be repeated on a finer grid in the 
vicinity of the found minimum. 

• * ■ 

Having found parameters a m i n , P m j n and y m m for the ascending branch we 

20 generate a dense t-grid {t s , s=l,2,*", N} and calculate the corresponding values of the 
interval duration as 

m * 

Ts-T^fyjl + Kia^fi^ia^&^t,-*^]). ' (12.26) 



25 A similar dense representation of the descending branch can be calculated in exactly 
the same way. The resulting bird-like curves are illustrated in panels A and C of 
Figure 20. The actual absolute and relative errors are indicated in the captions. The 
right hand side panels B and D represent hysteresis curves in two different 
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representations and each resulting from the curves shown in Panels A and C. The 
following computations of the hysteresis loop and its measure given by Eq. (10.17) is 
then performed as described in Example 10. 



5 EXAMPLE 13 

» • 

Creation of an RR-Hysteresis Loop 
With the Procedures of Examples 10-12 

In addition to the procedures generating a hysteresis loop on .the plane (QT- 
interval versus RR-interval) or an equivalent plane, one can introduce and assess a 
10 separate hysteresis of the duration of RR-interval versus exercise workload, which 
gradually varies, there-and-back, during the ascending and descending exercise 
stages. The RR-hysteresis can be displayed as loops on different planes based on just 

a single RR data set ,7^ } analysis. For example, such a loop can be displayed on 

the (X 1 , T^^plane, where T 1 = \t M — and is the time instant corresponding 
15 to the peak of the exercise load, or the center of the maximum load period, which may 

• « 

be determined according to numerical techniques described in the examples 10-12. 
The RR loop can also be introduced on the (W (t'^) 9 T^) or (T 1 , (T^)" 1 ) planes. 

Here JF(^)is a workload that varies . versus exercise stages, i.e. time and 

■ 

(T^ )' 1 represents the heart rate. 

20 In order to apply a numerical technique from the Example 1 0 (or any other as 

*■ • 

in Examples 1 1 or 12), one repeats essentially the same computational steps described 
in these examples. However, instead of considering both QT- {^ r >7^ r } and RR- 

{/^,T^} interval data sets, only a single RR data set is numerically processed 

> 

throughout the whole sequence of the described stages for a creation of a hysteresis 
25 loop. In this case variables t', (T^)' 1 or J^(?^)play the role of the second 

* 

alternating component that along with the first variable forms the RR - hysteresis 
plane. 
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The foregoing examples are illustrative of the present invention and are not to 
be construed as limiting thereof. The invention is defined by the following claims, 

• ■ 

with equivalents of the claims to be included therein. 
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THAT WHICH IS CLAIMED IS: 

1. A method of assessing cardiac ischemia in a subject to provide a measure 
. of cardiac or cardiovascular health in that subject, said method comprising the steps 

of: ' 

♦ 

5 (a) collecting a first QT- and RR- interval data set from said subject during a 

stage of gradually increasing heart rate; 

4 * 

(]b) collecting a second QT- and RR- interval data set from said subject during 

*r - 

1 

a stage of gradually decreasing heart rate, with said first and second QT- and RR- 
interval data sets being collected while minimizing the influence of rapid transients 
0 due to autonomic nervous system and hormonal influence on said data sets; 

(c) comparing said first QT- and RR- interval data set to said second QT- and 
RR- interval data set to determine the difference between said data sets; and 

(d) generating from said comparison of step (c) a measure of cardiac ischemia 
during stimulation in said subject, wherein a greater difference between said first and 

5 second data sets indicates greater cardiac ischemia and lesser cardiac, or 
cardiovascular health in said subject. 

i * 

2. The method according to claim 1, wherein said first and second QT- and 

r 

RR-interval data sets are collected under quasi-stationary conditions. 

0 

■ 

3. The method according to claim 1, wherein said stage of gradually increasing 
heart rate and said stage of gradually decreasing heart rate are each at least 3 minutes 
induration. 

5 4. The method according to claim 1, wherein said stage of gradually increasing 

heart rate and said stage of gradually decreasing heart rate are together carried out for 

» 

a total time of from 6 minutes to 40 minutes. 

5. The method according to claim 1, wherein: 

3 both said stage of gradually increasing heart rate and said stage of gradually 

■ 

decreasing heart rate are carried out between a peak rate and a minimum rate; and 
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a " 

» 

said peak rates of both said stage of gradually increasing heart rate and said 

< 

. stage of gradually decreasing heart rate are the same. . 

■ 

• ■ 

a 

6. The method according to claim 5, wherein: 

5 said minimum rates of both said stage of gradually increasing heart rate and * 

said stage of gradually decreasing heart rate are substantially the same. 

i * 

• » 

7. The method according to claim 1, wherein said stage of gradually 
decreasing heart rate is carried out at at least three different heart-rate stimulation 

■ 

10 levels. . ' • 

t 

■ 

8. The method according to claim 7, wherein said stage of gradually 

» 

increasing heart rate is carried out at at least three different heart-rate stimulation 

■ 

■ 

levels. - 

15 

> 

■ * 

9. The method according to claim 1, wherein said stage of gradually increasing 
heart rate and said stage of gradually decreasing heart rate are carried out sequentially 
in time. 

.20 10. The method according to claim 1, wherein said stage of gradually 

increasing heart rate and said stage of gradually decreasing heart rate are carried out 
separately in time. 

* 

11. The method according to claim 1, wherein said generating step is carried 
25 out by generating curves from each of said data sets. 

« 

a 

12. The method according to claim 1 1, wherein said generating step is carried 
out by comparing the shapes of said curves from data sets. ' 

. a 
■ 

* a 
• a 

30 13. The method according to claim 1 1, wherein said generating step is carried 

out by determining a measure of the domain between said curves. 
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14. The method according to claim 11, wherein said generating step is carried 
out by both comparing the shapes of said curves from data sets and determining a 
measure of the domain between said curves. 

* 

5 15. . The method according to claim 11, further comprising the step of 

■ 

* • • 

displaying said curves. 

■ 

■ 

16. The method according to claim 1, wherein said heart rate during said stage 
of gradually increasing heart rate does not exceed more than 120 beats per minute. 

t 

17. The method according to claim 1 , wherein said heart rate during said stage 
of gradually increasing heart rate exceeds 120 beats per minute. 

T 

1 8. The method according to claim 1, further comprising the step of: 

15 (e) comparing said measure of cardiac ischemia during stimulation to at least 

* a 

one reference value; and then - 

. (f) generating from said comparison of step (e) a quantitative indicium of 

cardiac or cardiovascular health for said subject. 

■ 

20 19. The method according to claim 18, further comprising the steps of: 

(g) treating said subject with a cardiovascular therapy; and then 
. (h) repeating steps (a) through (j) to assess the efficacy of said cardiovascular 
therapy, in which a decrease in the difference between said data sets from before said 
therapy to after said therapy indicates an improvement in cardiac health, in said 
25 subject from said cardiovascular therapy. 

• - 

20. The method according to claim 19, wherein said cardiovascular therapy is 
selected from the grojip consisting of aerobic exercise, muscle strength building, 
change in diet, . nutritional supplement, weight loss, stress reduction, smoking 
30 cessation, pharmaceutical treatment, surgical treatment, and combinations thereof. 
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21. The method according to claim 18, further comprising the step of 

■ 

assessing from said quantitative indicum the likelihood that said subject is at risk to 
experience a future ischemia-related cardiac incident. 

i 

5 22. A method of assessing cardiac ischemia in a subject to provide a measure 

of cardiac or cardiovascular health in that subject, said method comprising the steps 
of: 

(a) collecting a first QT- and RR- interval data set from said subject during a 
stage of gradually increasing exercise load and gradually increasing heart rate; and 

10 then, without an intervening rest stage, 

(b) collecting a second QT- and RR- interval data set from said subject during 
a stage of gradually decreasing exercise load and gradually decreasing.heart rate; 

(c) comparing said first QT- and RR- interval data set to said second QT- and 
RR- interval data set to determine the difference between said data sets; and 

15 (d) generating from said comparison of step (c) a measure of cardiac ischemia 

during exercise in said subject, wherein a greater difference between said first and 
second data sets indicates greater cardiac ischemia and lesser cardiac or 
cardiovascular health in said subject. 

• i 

• m 

20 23. The method according to claim 22, wherein said stage of gradually 

increasing exercise load and gradually decreasing exercise load are each at least 3 
minutes in duration. 

24. The method according' to claim 22, wherein said stage of gradually 
25 increasing exercise load and said stage of gradually decreasing exercise load are 

* 

* 

together carried out for a total time of from 6 minutes to 40 minutes. . 

« 

■ 

H • 

25. The method according to claim 22, wherein: 

both said stage of gradually increasing exercise load and said stage of 

• > 

30 gradually decreasing exercise load are carried out between a peak load and a 
minimum load; and 
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said peak loads of both said stage of gradually increasing exercise load and 
said stage of gradually decreasing exercise load are essentially the same. 

26. The method according to claim 25, wherein: 
5 said minimum loads of both said stage of gradually increasing exercise load, 

and' said stage of gradually decreasing exercise load are essentially the same. 



10 



27. The method according to claim 22, wherein said stage of gradually 
decreasing exercise load is canied out at at least three different load levels. 

28. The method according to claim 22, wherein said stage of gradually 
increasing exercise load is carried out at at least three different load levels. 



29. The method according to claim 22, wherein said stage of gradually 
15 increasing exercise load and said stage of gradually decreasing exercise load are 
carried out sequentially in time. 



30. The method according to claim 22, wherein said stage of gradually 
increasing exercise load and said stage of gradually decreasing exercise load are 
20 carried out separately in time. 



• 31 . The method according to claim 22, wherein said generating step is carried 
out by generating curves from each of said data sets. 



25 32. The method according to claim 31, wherein said generating step is carried 

out by comparing the shapes of said curves from said data sets. 



33. The method according to claim 31, wherein said generating step is carried 
out by determining a measure of the domain between said curves. 

30 



* 
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* 

34. The method according to claim 31, wherein said generating step is carried 
out by both comparing the shapes of said curves from data sets and determining a 
measure of the domain between said curves. 

5 .35. The method according to claim 31, further comprising the step of 

displaying said curves. 

- 

36. The method according to claim 22, wherein said* heart rate during said 
stage of gradually increasing exercise load does not exceed more than 120 beats per 

10 minute. 

9 

37. The method . according to claim 22 r wherein said heart rate during said 
stage of gradually increasing heart rate exceeds 120 beats per minute. 

* 

* ■ 

■ 

■ 

15 38. The method according to claim 22, further comprising the step of: 

(e) comparing said measure of cardiac ischemia during exercise to at least one 
reference value; and then 

■ 

(f) generating from said comparison of step (e) a quantitative indicium of 
cardiac or cardiovascular health for said subject. 

* 

20 

39. The method according to claim 38, further comprising the steps of: 

(g) treating said subject with a cardiovascular therapy; and then 

* 

(h) repeating steps (a) through (f) to assess the efficacy of said cardiovascular 
therapy, in which a decrease in the difference between said data sets from before said 

25 therapy to after said therapy indicates an improvement in cardiac or cardiovascular 
health in said subject from said cardiovascular therapy. 

■ 

4 

40. The method according to claim 39, wherein said cardiovascular therapy is 
selected from the group consisting of aerobic exercise, muscle strength building, 

30 change in diet, nutritional supplement, weight loss, stress reduction, smoking 
cessation, pharmaceutical treatment, surgical treatment, and combinations thereof. 
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» 

41. The method according to claim 38, further comprising the step of 
assessing from said quantitative indicium the likelihood that said subject is at risk to 
experience a future ischemia-related cardiac incident. 

* ■ 

■ * 

.5 • 42. A method of assessing cardiac ischemia in a subject to provide a measure 

of cardiac or cardiovascular health in that subject, said method comprising the steps 
of: 

(a) collecting a first QT- and RR- interval data set from said subject during a 
stage of gradually increasing exercise load and gradually increasing heart rate at least 

10 5 minutes in duration, and at at least three different load levels; and then, without an 

■ 

intervening rest stage; 

(b) collecting a second QT- and RR- interval data set from said subject during 
a stage of gradually decreasing exercise load and gradually decreasing heart rate at 
least 5 minutes in duration, and at at least three different load levels; 

15 (c) comparing said first QT- and RR- interval data set to said second QT- and 

RR- interval data set to determine the difference between said data sets; 

* X 

(d) generating from said comparison of step (c) a measure of cardiac ischemia 

» 

during exercise in said subject, wherein a greater difference between said first and 
second data sets indicates greater cardiac ischemia and lesser cardiac or 
20 cardiovascular health in said subject; 

(e) comparing said measure of cardiac ischemia during exercise to at least one 
reference value; and then 

(f) generating from said comparison of step (e) a quantitative indicium of 
cardiac or cardiovascular health for said subject. 

25 . 

* 

43. The method according to claim 42, wherein said stage of gradually. 

■ 

increasing exercise load and said stage of gradually decreasing exercise load are 
together carried out for a total time of from 6 minutes to 40 minutes. 

* 

30 44. The method according to claim 42, wherein; 
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both said stage of gradually increasing exercise load and said stage of 
gradually decreasing exercise load are carried out between a peak load and a 
minimum load; and 

said peak loads of both said stage of gradually increasing exercise load and 
5 said stage of gradually decreasing exercise load are essentially the same. 

, * 

• * 

45. The method according to claim 44, wherein: 

said minimum loads of both said stage of gradually increasing exercise load 
and said stage of gradually decreasing exercise load are essentially the same. 

10 

46. The method according to claim 42, wherein said stage of gradually 
increasing exercise load and said stage of gradually decreasing exercise load are 
carried out sequentially and contiguously in time. 

■ 

15 47. The method according to claim 42, wherein said stage of gradually 

increasing exercise load and said stage of gradually decreasing exercise load are 
carried out separately in time. 

» * 

48. The method according to claim 42, wherein said generating step is carried 
20 out by generating curves from each of said data sets. 

■ 

* 

49. The method according to claim 48, wherein said generating step is carried 

■ 

* 

out by comparing the shapes of said curves from said data sets. 

25 50. The method according to claim 48, wherein said generating step is carried 

out by determining a measure of the domain between said curves. 

■ 

* ■ 
* 

r • i 

' 51 . The method according to claim 48, wherein said generating step is carried 
out by both comparing the shapes of said curves from said data sets and determining a m 
30 ■ measure of the domain between said curves. 
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' 52. The method according to claim 48, further comprising the step of 
displaying said curves. 

- 

53. The method according to claim 42, wherein said heart rate during said 
5 stage of gradually increasing exercise load does not exceed more than 120 beats per 

minute. 

54. The method according to claim 42, wherein said heart rate during said 
stage of gradually increasing heart rate exceeds 120 beats per minute. 

10 

* 

55. The method according to claim 42, further comprising the steps of: 

(g) treating said subj ect with a cardiovascular therapy; and then 

(h) repeating, steps (a) through (f) to assess the efficacy of said cardiovascular 
therapy, in which a decrease in the difference between said data sets from before said 

15 therapy to after said therapy indicates an improvement in cardiac or cardiovascular 

* 

health in said subject from said cardiovascular therapy. 

» ♦ 

■ 

56. The method according to claim 55, wherein said cardiovascular therapy is 
selected from the group consisting of aerobic exercise, muscle strength building, 

20 change in diet, nutritional supplement, weight loss, stress reduction, smoking 
cessation, pharmaceutical treatment, surgical treatment, and combinations thereof. 

* 

m 

t i 

- m 

57. The method according to claim 42, further comprising the step of 
assessing from said quantitative indicium the likelihood that said subject is at risk to 

25 experience a future ischemia-related cardiac incident. 

* 

• • • 

58. A method of assessing cardiac ischemia in a subject to provide a measure 
of cardiac or cardiovascular health in that subject, said method comprising the steps 

* 

of: 

30 (a) collecting a first QT- and RR- interval data set from said subject during a 

• • * 

stage of gradually increasing heart rate, said heart rate being gradually increased in 
response to actual heart rate data collected from concurrent monitoring of said patient; 
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I 

* • 

(b) collecting a second QT- and RR- interval data set from said subject during 

♦ 

a stage of gradually decreasing heart rate, said heart rate being gradually decreased in 
response to actual heart rate data collected from concurrent monitoring of said patient; 

(c) comparing said first QT- and RR- interval data set to said second QT- and 

■ 

RR- interval data set to determine the difference between said data sets; and 

(d) generating from said comparison of step (c) a measure of cardiac ischemia 

* 

during stimulation in said subject, wherein a greater difference between said first and 
second data sets indicates greater, cardiac ischemia and lesser cardiac or 
cardiovascular health in said subject. ' 



.59. The method according to claim 58, wherein said first and second QT- and 
RR-interval data sets are collected without an intervening rest period. 

60. The method according to claim 58, wherein said stage of gradually 
1 5 increasing heart rate and gradually decreasing heart rate are each at least 3 minutes in 

* 

duration. 



61. The method according to claim 58, wherein said stage of gradually 

• • • 

increasing heart rate and said stage of gradually decreasing heart rate are together 
20 carried out for a total time of from 6 minutes to 40 minutes. . 



62. The method according to claim 58, wherein: 
' both said stage of gradually increasing heart rate and said stage of gradually 
decreasing heart rate are carried out between a peak rate and a minimum rate; and 

* 

25 said peak rates of both said stage of gradually increasing heart rate and said 

* 

■ 

stage of gradually decreasing heart rate are substantially the same. 



63. The method according to claim 62, wherein: 

» * 

said minimum rates of both said stage of gradually increasing heart rate and 
30 said stage of gradually decreasing heart rate are the same. 
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64. The method according to claiin 58, wherein said stage of gradually 

4 * 

decreasing heart rate is carried out at at least three different heart-rate stimulation 
levels. 



5 65. The method according to claim 64, wherein said stage of gradually 

increasing heart rate is carried out at at least three different heart-rate stimulation 
levels. 



66. The method according to claim 58, wherein said stage of gradually 
1 0 increasing heart rate and said stage of gradually decreasing heart rate are carried out 

sequentially in time. 

■ 

67. The method according to claim 58, wherein said stage of gradually 
increasing heart rate and said stage of gradually decreasing heart rate are carried out 

1 5 separately in time. 



68. The method according to claim 58, wherein said generating step is carried 
out by generating curves from each of said data sets. 



20 69. The method according to claim 68, wherein said generating step is carried 

r 

out by comparing the shapes of said curves from said data sets. 

■ i 

70. The method according to claim 68, wherein said generating step is carried 
out by determining a measure of the domain between said curves. 



71. The method according to claim 68, wherein said generating step is carried 

■ 

out by both comparing the shapes of said curves from said data sets and detennining a 
measure of the domain between said curves. 



* r 

30 72. The method according to claim 68, further comprising the step of 

displaying said curves. 
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73. The method according to claim 58, wherein said heart rate during said 
stage of gradually increasing heart rate does not exceed more than 120 beats per 
'minute. 

\ 

* 

* 

* • > 

5 74. The method according to claim 58, wherein said heart rate during said 

stage of gradually increasing heart rate exceeds 1 20 beats per minute. 

■ 

75. The method according to claim 58, further comprising the step of: 

(e) comparing said measure of cardiac ischemia during stimulation to at least 
10 one reference value; and then 

(f) generating from said comparison of step (e) a quantitative indicium of 
cardiac or cardiovascular health for said subject. 

* ■ 

« 

76. The method according to claim 75, further comprising the steps of: 
15 . (g) treating said subject with a cardiovascular therapy; and then 

■ 

(h) repeating steps (a) through (f) to assess the efficacy of said cardiovascular 
therapy, in which a decrease in the difference between said data sets from before said 
therapy to after said therapy indicates an improvement in cardiac or cardiovascular 
health in said subject from said cardiovascular therapy. 

20 

■ 

77. The method according to claim 76, wherein said cardiovascular therapy is 
selected from the group consisting of aerobic exercise, muscle strength building, 
change in diet, nutritional * supplement, weight loss, stress reduction, smoking 

* 

cessation, pharmaceutical treatment, surgical treatment, and combinations thereof. 

> 

25 

. 78. The method according to claim 75,- further, comprising the step of 
assessing from said quantitative indicium the likelihood that said subject is at risk to 
experience a future ischemia-related cardiac incident. 

» 

30 79. A method of assessing cardiac ischemia- in a subject to provide a measure 

♦ 

of cardiac or cardiovascular health in that subject, said method comprising the steps, 
performed on a computer system, of: 
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(a) providing a first QT- and RR- interval data set collected from said subject 

> ■ 

during a stage of gradually increasing heart rate; 

(b) providing a second QT- and RR- interval data set collected from said 
subject during a stage of gradually decreasing heart rate; 

■ 

5 (c) comparing said first QT- and RR- interval data set to said second QT- and 

RR- interval data set to determine the difference between said data sets; and 

(d) generating from said comparison of step fcja measure of cardiac ischemia 
during stimulation in said subject, wherein a greater difference between said first and 
second data sets indicates greater cardiac ischemia and lesser cardiac or 
10 cardiovascular health in said subject. 

80. The method according to claim 79, wherein said first and second QT- and 
RR-interval data sets are collected while minimizing the influence of rapid transients 
due to autonomic nervous system and hormonal control on said data sets 

« * 

81. The method according to claim 79, wherein said first and second QT- and 
RR-interval data sets are collected without an intervening rest stage. 

82. The method according to claim 79, wherein said generating step is carried 

* 

20 out by generating curves from each of said data sets. 

■ 

83. The method according to claim 79, wherein said generating step is carried 
out by comparing the shapes of said curves from data sets. 

■ 

» 

25 > 84. The method according to claim 79, wherein said generating step is carried 

out by determining a measure of the domain between said curves. 

■ 

85. The method according to claim 79, wherein said generating step is carried 
out by both comparing the shapes of said curves from data sets and determining a 
30 measure of the domain between said curves. 
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86. The method according to claim 79, further comprising the step of 
displaying said curves. 



87. The method according to claim 79, wherein said comparing step is carried 
5 . out by: , 

(i) filtering said first and second QT- and RR- interval data sets; 

(ii) generating a smoothed hysteresis loop from said filtered first and second 

i <- 

QT- and RR- interval data sets; and then 

(Hi) determining a measure of the domain inside said smoothed hysteresis 

10 loop. 



88. The method according to claim 79, wherein said comparing step is carried 

. out by: 

(i) filtering said first and second QT- and RR- interval data sets; 

* ■ 

15 (ii) generating preliminary minimal values for said first and second QT- and 

RR- interval data sets; 

(Hi) correcting said preliminary minimal values; 
. .(iv) generating first and second preliminary smoothed curves from each of said 
filtered data sets; 
20 (v/ correcting said preliminary smoothed curves; 

■ 

(vi) fitting said preliminary smoothed curves; 

(vii) generating a smoothed hysteresis loop from said first and second fitted 
smoothed curves; and then 

(viii) determining a measure of the domain inside said hysteresis loop. 

25 

89. The method according to claim 79, wherein said comparing step is carried 

out by: 

* 

(i) filtering said first and second QT- and RR- interval data sets by moving 

•> 

average smoothing; 

*~ . 

30 ■ ' (ii) generating a smoothed hysteresis loop from said filtered first and second 
QT- and RR- interval data sets; and then 

* 

(Hi) determining a measure of the domain inside said hysteresis loop. 
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90. The method according to claim 79, wherein said comparing step is carried 

out by: 

5 (i) sequentially combining said first and second QT- interval data sets into one 

■ 

■ 

combined QT-interval set; 

fii) sequentially combining said first and second RR- interval data sets into 
one combined RR-interval set; 

(Hi) filtering said combined QT- and RR- interval data sets; 
10 (iv) evaluating coordinates of the minima of said combined and filtered QT- 

and RR- interval data sets; 

(v) defining central data subsets including the minima of said combined and . 
filtered QT- and RR-interval data sets; 

(vi) fitting each said central data subset with a parabola that has a minimum 
1 5 inside said data subset; 

(vii) determining coordinates of the minimum for each said parabola; 

(viii) defining reduced data subsets by excluding said central data sets from 
said combined and filtered QT- and RR-interval data sets; 

fix) generating first transformed QT- and RR -interval data sets by mapping 
20 falling off branches of said reduced data sets by a standard monotonously decreasing 
mapping function of RR- or QT-intervals, which depends on at least one discrete or 
continuous parameter; 

(x) generating second transformed QT- and RR -interval data sets by mapping 
growing branches of said reduced data sets by a standard monotonously increasing 
25 mapping function of RR- or QT-intervals, which depends on at least one discrete or 

■ 

continuous parameter; 

fxi) fitting said first and second transformed sets by a standard function of 
time depending on one or more parameters whose values are determined by the ' 

* 

condition of the best fit; 
30 fxii), completing said fitting step by determining the parameters of said 

mapping functions of RR- and QT-intervals by the condition of the best fit for each 
branch of said QT- and RR -interval data sets; 
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(xiii) generating a hysteresis loop from said best fit functions; and then 

(xiv) evaluating a measure of the domain inside said hysteresis loop. 

* 

91. A computer system for assessing cardiac ischemia in a subject to provide a 
5 measure of cardiac or cardiovascular health in that subject, said system comprising: 

(a) means for providing a first QT- and RR- interval data set from said subject 

■ 

during a stage of gradually increasing heart rate; 

(b) means for providing a second QT- and RR- interval data set from said 
subject during a stage of gradually decreasing heart rate; 

* • • 

10 • (c) means for comparing said first QT- and RR- interval data set to said second 

» 

QT- and RR- interval data set to determine the difference between said data sets; and 

(d), means for generating from said comparison of step (c) a measure of 
cardiac ischemia during stimulation in said subject, wherein a greater difference 
between said first and second data sets indicates greater cardiac ischemia and lesser 
15 cardiac or cardiovascular health in said subject. 

* 

92. The system according to claim 91, wherein said generating is carried out 
by generating curves from each of said data sets. 

20 93. The system according to claim 91, wherein said generating is carried out 

by comparing the shapes of said curves from data sets, 

94. The system according to claim 91, wherein said generating is carried out 
by determining a measure of the domain between said curves. 

25 

95. The system according to claim 91, wherein said generating is carried out 
by both comparing the shapes of said curves from data sets and determining a 

* * « ♦ 

measure of the domain between said curves. 

4 

* 

* 

• * 

30 96. A computer program product for assessing cardiac'ischemia in a subject to 

provide a measure of cardiac or cardiovascular health in that subject, said' computer 
program product comprising a computer usable storage medium having computer 
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readable program code means embodied in the medium, the computer readable 
program code means comprising: 

- (a) computer readable program code means for comparing a first QT- and RR- 
interval data set to a second QT- and RR- interval data set to determine the difference 
5 between said data sets; and 

(b) computer readable program code means for generating from said 

* 

comparing of step (b) a measure of cardiac ischemia during stimulation in said 
subject, wherein a greater difference between said first and second data sets indicates 
greater cardiac ischemia and lesser cardiac or cardiovascular health in said subject. 

10 

97. The product according to claim 96, wherein said generating is carried out 
by generating curves from each of said data sets. 

98. The product according to claim 96, wherein said generating is carried out 

« 

15 by comparing the shapes of said curves from data sets. 

♦ 

99. The product according to claim 96, wherein said generating is carried out 

4 

by determining a measure of the domain between said curves. 

n 

20 100. The product according to claim 96, wherein said generating is carried out 

by both comparing the shapes of said curves from data sets and determining a 
measure of the domain between said curves. 
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